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Abstract 



The transverse momentum of a colour- singlet massive particle in a hadronic 
5^ \ collision is built up by two components, the perturbative effect of parton 

branchings and the nonperturbative effect of primordial k±. In previous 
studies of transverse momentum spectra for Z° production at the Tevatron, 
the best fit to the experimental data are given when the primordial k± is set 
to a much higher value than what is expected considering the confinement 
of partons in the proton. We here investigate the possibility that the reason 
for this is that too few branchings are generated in showers, compared to 
the evolution used in the tuning of parton densities. This could then be 
compensated by increasing the value of Aqcd- The study is done using 
the regular Pythia showering routines and a new algorithm where the 
branchings are ordered in transverse momentum instead of virtuality Q^. 



1 Introduction 



Again, of bodies some are composite, others the elements of which these 
composite bodies are made up. These elements are indivisible and unchange- 
able, and necessarily so, if things are not all to be destroyed and pass into 
non-existence, but are to be strong enough to endure when the composite bod- 
ies are broken up, because they possess a solid nature and are incapable of 
being anywhere or anyhow dissolved. It follows that the first beginnings must 
be indivisible, corporeal entities. 

Epicurus "Letter to Herodotus" 
approximately 300 B.C. 

What are the building blocks of the universe? That is a very profound question, 
pondered upon by philosophers through the milleniums. For a long time the theories were 
purely metaphysical, with no possibility whatsoever to be tested experimentally. Thanks 
to major breakthroughs in theoretical and experimental physics during the twentieth 
century we now have the Standard Model, describing those elementary particles and their 
interactions. This model has been extremely successful in predicting the outcome of 
experiments, performed at increasingly large accelerators. 

The elementary particles are divided into quarks, leptons and gauge bosons. Also the 
corresponding antiparticles exist. The gauge bosons mediate the four fundamental forces 
of nature: electromagnetism, the weak force, the strong force and gravity. The interesting 
particles for this thesis are the quarks and the gluons, the gauge bosons mediating the 
strong force. These are denoted q and g, with q denoting an antiquark. The theory 
describing the interactions of quarks and gluons is called quantum chromo dynamics, 
abbreviated QCD. 

Free quarks have never been observed. The theoretical explanation for this has to 
do with the behaviour of the strong force. Here, it is useful to make a comparison with 
electromagnetism. Instead of electromagnetic charge, the strong force acts on particles 
with colour charge. This is, of course, not real colour, but is named this way since the 
behaviour of the colour charges resembles the behaviour of real colours in the way that 
mixing all the colour charges makes a colour neutral object. Of the elementary particles 
it is the quarks and gluons that carry colour charge. The crucial point is that the force 
between two colour charged particles does not decrease with the distance separating them, 
as is the case for electromagnetism. This means that it takes an infinite amount of energy 
to separate a qq pair, if they were to be connected to each other by a colour field. When 
the distance between the quarks gets sufficiently large, the potential energy becomes large 
enough to create a new qq pair. The new qq pair screens the colour field between the 
original charges. But now there arc two colour connected qq pairs. So when the quark 
and the antiquark of a colour connected qq pair are moved away from each other, new qq 
pairs are created, but free quarks will not be obtained. 

Another important difference from electromagnetism is that the particle mediating the 
strong force, the gluon, is not a colour neutral object. This means that a gluon can emit 
a gluon. So there are three different QCD branchings of particles: q — > qg and g — > qq, 
which have their corresponding branchings in electromagnetism (q and 7 — > qq 

respectively), and g — > gg, which has no similar branching in electromagnetism. 

Quarks can carry the colour charges red, green or blue, the antiquarks the complemen- 
tary anticolours and gluons one colour and one anticolour. As noted above only colour 
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neutral objects can appear as free particles. Colour neutral particles out of quarks can be 
obtained in two different ways, either qq (colour-anticolour) or qqq (the mixture of red, 
green and blue makes a colour neutral object). The former particles are called mesons 
and the latter baryons (of course also antibaryons qqq exist), of which the proton and 
the neutron are the most well-known. The common name for all particles consisting of 
quarks is hadrons. 



quarks are fairly light and will be assumed to be massless in this thesis, the c and b 
quarks are considerably heavier and the t quark is very heavy and too shortlived to be 
bound in any hadron. 

To test the Standard Model experimentally, particles are accelerated to high energies 
and then collided. In the collisions the energy can be used to create new particles. In this 
way many particles predicted by the Standard Model have been found. Much of the effort 
in particle physics today is aimed at finding the Higgs particle, the last particle predicted 
by the Standard Model that has not yet been found. This search has been performed 
at the Tevatron, an accelerator at Fermilab outside of Chicago, and will be continued at 
LHC (Large Hadron Collider), an accelerator under construction at CERN in Geneva. 

At the Tevatron protons and antiprotons are collided and at LHC protons are collided 
with protons. The energies are high enough to resolve the quarks and the gluons in the 
proton; these are called partons with a common name. The process of interest is often 
a hard collision between partons, which is not easy to extract from the experimental 
data or calculate theoretically. Firstly, in the final state, the partons created at the 
high energy scale of the hard collision branch into showers of partons. These parton 
showers evolve to a lower energy scale, where the partons are bound into hadrons. It is 
these hadrons that show up in the detector. The hadronization of partons can not be 
described by perturbative methods. Secondly, in the initial state, when the accelerated 
particles are protons, there is a problem of knowing the inner structure of the proton at 
different energy scales. High energy scales correspond to short distance interactions, as 
hard collisions of partons. Since the process of interest is a hard collision of partons, the 
distribution of partons in the proton at this energy scale must be known. The partons of 
the hard coUision can be generated by letting the partons of the proton at a lower energy 
scale develop initial-state parton showers up to the scale of the hard collision. These 
initial-state parton showers will be the subject of the study in this thesis. 

At the scale of particle physics our universe is described by quantum mechanics, which 
means that it is not possible to predict the outcome of a single event. Only the probabili- 
ties of what is going to happen can be predicted. It is therefore necessary to average over 
a large number of events in a meaningful way. Perturbative methods can be applied to 
describe events down to a low energy scale, under which nonperturbative methods have 
to be used. In an event hundreds of particles can be produced. This means that even if 
the calculations could be done using perturbative methods, the complexity would exclude 
many approaches. One approach that does work is doing Monte Carlo simulations of a 
large number of events. Both initial-state and final-state parton showers are generated. 
This way spectra of different variables, that can be compared to experiment, are obtained. 
Event generators are used to interpret the outcome of experiments and to predict what 
to look for in the future. One of the most successful event generators, that will be used 
in this thesis, is called Pythia [1, 2]. 

One interesting variable to simulate in event generators is transverse momenta p±, the 
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momentum that particles have transverse to the accelerated beam. In the branchings the 
partons undergo, they get recoil transverse to the beam axis. This means that the more 
branchings, the more p± of the partons. However there is also an opposite effect, the 
accumulated of a parton is split between the two new partons of a branching. Still, 
the branchings affect pj_, so to look at p± is a good way to find out something about the 
branchings of an event. p± can be seen as a measure of how violent a process is. 

One process where it is possible to make a comparison between theory, event generators 
and experimental data in a rather clean way is qq — > Z°. When the Z° decays to an e+e" 
pair, its p± is reasonably easy to measure. The p± of the Z'^ boson is the sum of the 
p± of the partons in the hard collision, where it was created. The p± of these partons 
is built up by two components; a perturbative component from the parton branchings 
and a non-perturbative component called primordial k± to be described as follows. The 
partons in the proton have an uncertainty in momenta that, according to Heisenberg's 
uncertainty relation, should be inversely proportional to the radius of the proton. This 
gives rise to a transverse momentum, that is called primordial k±. For high p± values of 
the Z° boson the perturbative component from the parton branchings dominates, while 
the low p± part of the p^-spectrum is sensitive to primordial k±. 

The p± of the Z*^ and bosons in the process qq — ^ Z^/W*^ have been measured 
at the Tevatron [3, 4, 5] and there is a problem with the Pythia prediction of the p±- 
spectrum for small p± as follows. To get the best fit to experimental data for small p± the 
value of the primordial k± has to be set considerably higher than expected, considering 
the radius of the proton. 

This thesis investigates one possible cause to this problem. The study will be made 
for the case of Z° production. In Pythia, some perturbative branchings, that normally 
are assumed to contribute, are not simulated. This means that the p_L-spectrum of Z° 
boson will be shifted to lower values of p±. To compensate for this the other component 
building up the p± of the Z° boson, the primordial k±, has to be increased. 

If the cause of the problem is that Pythia generates too few branchings, there might 
be a way of solving it. In the expression for the running strong coupling constant ctsiQ^) 
there is a parameter Aqcd that has to be measured experimentally. An increase of Aqcd 
increases the rate of branchings. At first sight this does not seem to provide any solution 
to the problem, since the same Aqcd used to determine the parton densities, which are 
used in Pythia, must be used in the Pythia parton showers. However, the parton 
densities are leading log. In Pythia there arc some corrections to pure leading log 
DGLAP evolution, introduced to give better agreement with next to leading log results. 
It is these corrections that decrease the rate of branchings, so to compensate for this Aqcd 
can be increased. What is studied in this thesis is how much Aqcd should be increased to 
compensate for the corrections and to what extent this solves the problem with primordial 
k^. 

In Section 2 the concepts of parton showers and the DGLAP equations are introduced 
and a more specific description of the way of doing initial-state parton showers in Pythia 
is given. Section 3 is a description of the problem with the too high value of primordial 
kj_, when comparing to experiment. The corrections to DGLAP evolution in Pythia are 
described in Section 4. The study of the generation of p_|_-spectra is presented in Section 5. 
In this section also a new algorithm for doing initial-state showers is introduced. This 
new algorithm does not include heavy quarks. Section 6 is a attempt to find a way of 
introducing heavy quarks in the new algorithm. And finally Section 7 is a summary of 
the results of this thesis. 
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Figure 1: From the three processes q —>■ qg, g ^ gg and g — * qq it is possible to build up 
complicated showers. The emissions are ordered in Q^. 

2 Parton showers 

To simulate events, that can be compared with the data from accelerators, either matrix 
element calculations or parton showers can be used. In the matrix element approach 
amplitudes and phases are used to do a fully quantum mechanical calculation. Calcu- 
lating matrix elements order by order in the strong coupling constant ag is a good, but 
complicated, method to get cross sections for processes. As each branching contributes 
with one order in as, the matrix element calculations get too complicated for processes 
involving more than a few branchings. An approximate approach is needed. In such an 
approach probability distributions, instead of amplitudes and phases, can be used. This 
is an ideal method for doing computer simulations. Using the probabilities for the differ- 
ent branchings, given by the DGLAP equations [6], one parton branching into a shower 
of partons can be simulated [7]. This method can be applied to events with arbitrary 
many branchings, but is only accurate to what is called leading log, where to each order 
in as, only the most divergent contribution is described correctly. In this parton shower 
approach the histories of the partons are given, so they can be traced back through the 
branchings. 

2.1 General concepts of parton showers 

In the introduction the proton is said to consist of three quarks; this is a simplification. 
In the proton different processes is constantly occuring: q — > qg, g — gg, g — > qq and 
their inverses. The quantum numbers of a quark originating from a gluon splitting are 
always compensated by the quantum numbers of the antiquark originating from the same 
splitting. These quarks are called sea quarks, while the three quarks that determine the 
quantum numbers of the proton are called valence quarks in analogy with atomic physics. 

As a consequence of these constantly occuring processes, the structure of the proton 
depends on the momentum scale it is probed by. A higher Q"^ can resolve more gluons 
and sea quarks. To describe the structure of the proton, parton density functions fa{x, Q^) 
are used. They express the probability of finding a parton of type a with a fraction x 
of the total proton momentum in a proton probed with Q^. To calculate parton density 
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functions nonperturbative input is needed. This nonperturbative input is equivalent to 
having the /a(x, Q"^) defined at one Q"^ — Qq. Given such initial conditions, the parton 
density function at another Q"^ can be calculated, see Subsection 2.2. The dependance 
of the parton density functions is often called scaling violations, since the /„ are no longer 
functions of x alone. 

Another way of describing this gradual resolution of the proton is to start with one 
initial parton and do an evolution in Q^. As increases the parton branches into more 
and more partons and a shower of partons arises. Such a parton shower is built up by 
the three processes (if processes involving photons are neglected) q qg, g — gg and 
g — > qq, see Fig. 1. To do a simulation of parton showers, the probability for these 
processes as a function of is needed. In Subsection 2.2 the DGLAP equations, giving 
these probabilities, are introduced. 

This parton shower approach can be used to generate events at particle colliders. An 
event is separated into one hard process, with the highest Q^, and many softer interactions 
both before and after the hard process. To describe these softer interactions, the parton 
shower approach is useful. The parton showers of an event can be separated into final-state 
showers and initial-state showers. The final-state showers evolve if the final particles of the 
collision are quarks or gluons. This can be simulated by an evolution in Q^, starting with 
the quarks at the scale of the hard collision and going down to the scale of hadronization. 
If the colliding particles are protons there is also initial-state parton showers, starting 
with the partons at some low scale Qq going up to the scale of the hard collision. In Fig. 1 
the parton shower starts with a quark 1 at Qo) which branches into a quark 2 and a gluon 
3 at a higher Q^. At an even higher quark 2 branches into a quark 4 and gluon 5 and 
at another gluon 3 splits into a quark 6 and an antiquark 7. The branchings continue 
up to the scale of the hard collision. Only one of the partons in the cascade takes part 
in the hard collision. When the partons that participates in the hard collisions are given, 
they can be traced back to the low scale by backwards evolution, to be described below. 
The backwards evolution starts with the parton of the hard collision, for instance gluon 
8, tracing it to the branching of gluon 5 at a lower Q^. Gluon 5 is emitted from quark 
2 at an even lower Q^. Quark 2 originates from the branching of quark 1, which is the 
initial parton of the shower at scale Qq. In this way the p_[_ of the colliding particles, and 
consequently the p±_ of the produced particle, can be calculated. 



2.2 DGLAP equations 

In DGLAP evolution the emissions are ordered in spacelike virtualities and the prob- 
ability of the process a — > 6c is given by 

dVa^hc = ^Pa-.bc{z) (1) 

where Pa^bc{z) is the splitting function depending on z, the fraction of the momentum of 
particle a taken by particle h. For the processes that are used to build up a QCD shower, 
the splitting functions are given by 



3 1-2 

z{l - z) 



(2) 
(3) 
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P.-^M - liz' + il-zf) (4) 



In eq. (1) the running strong coupling constant is used 



(33 - 2nf ) In (^) 

where nf is the number of flavours and Aqcd is a parameter of the Standard Model that 
has to be measured experimentally. Aqcd will be of great importance for this thesis. 

The parton shower approach is an exclusive approach where one parton, that branches 
into a shower as increases, is studied. To include all the partons of the proton, one can 
study the behaviour of parton densities as increases. Beginning with initial conditions 
for the parton densities at one Qq, considering the evolution of all the partons, the parton 
densities at another can be calculated. The dependence on of the parton density 
for a particle of type b at the value x is given by 

^^ = ^f^-^LM-')^P.^J^] (6) 

d(ln(52) ^Jx y 27r \y J 

2.3 Initial-state showers in Pythia 

In this thesis initial-state showers in Pythia [8, 9, 10] are studied. Pythia uses the 
parton shower approach to do Monte Carlo simulations of events. The DGLAP equations 
are used to simulate at which Q^, and with which z, the branchings occur. 

In the evolution a branching can only occur if the parton has not already branched. 
This introduces a Sudakov form factor to the expression for the probability 



To do a Monte Carlo simulation, a function N{Q'^) for the number of partons that has 
not branched is needed 



N{Q')^N{Ql){l- ll dV) (8) 

a random number R times the 
N{Ql). So the at which the branching occurs is given by 



'Ql 

This number is then set to a random number R times the number of partons at Qq, 



/ dV ^ 1-R^R 

using the fact that 1 — i? is a random number between and 1 if is it. The z of the 
branching can then be generated by setting 

/' dz'Pa^bc{z') ^rC dz'Pa-.bc{z') (10) 
JO JO 
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As seen in eq. (2) and eq. (3) the probabihty of emissions of very soft gluons, where the 
gluon takes a very small part of the momentum, goes to infinity. It is therefore necessary 
to introduce a cut-off e in z. In the case of q — > qg only an upper cut-off is needed, while 
in the case of g — gg both an upper and a lower cut-off is needed. Instead of integrating 
from on the left side and up to 1 on the right side in eq. (10), the integration is from e 
and up to 1 - e respectively. 

In the forward evolution is generated this way for the different possible branchings 
(q qg, g ^ gg and g — > qq) and the one with the lowest is chosen. Forward evolution 
can be very inefficient when a specific kind of hard process is studied, since the final-state 
partons are not known beforehand. If no final-state partons with the right properties 
to participate in the hard process are generated in the evolution, it can not be used. A 
more efficient method, used in Pythia, is to fix the hard scattering using evolved parton 
densities and thereafter use backwards evolution to form this into exclusive events. Such 
an approach can again be formulated in Monte Carlo terms, but now starting from eq. (6) 
instead of eq. (1). This gives 

where z has been rewritten as x, the momentum fraction of the parton before (in the 
ordering of the backwards evolution) the branching, divided by y, the momentum fraction 
of the parton after the branching. In the exponential the probability is integrated from 
Q'^ to Qmax) the virtuality of the preceding step, instead of from Ql to Q^, as is the case 
for forward evolution. 

The same procedure can then be repeated starting at the of the branching. In 
this way the partons of the hard interaction are evolved back to the initial partons at the 
lower scale Ql- 



3 Primordial k± 

The partons are confined in the proton of radius 0.7 fm. According to Heisenberg's 
uncertainty relation this means that they have a momentum inside the proton. The part 
of this momentum that is in the beam direction gives the fraction x of the total proton 
momentum carried by the parton, described by parton density functions. But there is 
also a part of this momentum that is transverse to the beam direction. This transverse 
momentum of the partons, due to the confinement in the proton, is called primordial k± 
and its mean value can be estimated by 

, h 0.2 GeV -fm _ , , 

(kj) ft; — «i — — ft 0.3 GeV 12 

Tp 0.7 fm 

To compare this value to experiment, we can look at the ]9^-spectrum of the produced 
Z° boson in the process qq — > Z°. The reason for looking at this particular process is 
that when the Z° boson decays to a e~^e~ pair the p± of the boson is reasonably easy 
to measure. The p^-spectrum of the Z° boson depends on the value of the primordial 
k± and on the p^-kicks the quarks get when they emit gluons. For high p± values of 
the Z° boson the effect of the p_L-kicks in the shower dominates, while the low p± part 
of the p_L-spectrum is sensitive to primordial k±. Fig. 2 gives a schematical picture of 
this. In one set of data from the Tevatron the Pythia prediction of the p_L-spectrum 
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Primordial k± 



Primordial k± 




Figure 2: Production of a Z^'-boson in a parton shower model. At each branching the 
accumulated p± is split between the particles in 2; vs. 1 — z fractions. New p± values 
for the partons after the branching are generated, and the shower continues. This means 
that showers with many branchings are not very sensitive to the value of primordial k±. 

fits the experimental data best for a primordial k± value of 2.15 Gev [11]. This value is 
considerably higher than the value estimated by Heisenberg's uncertainty relation. One 
reason for this could be that Pythia has too low a rate of emissions. Then there would 
be too few p^-kicks and the value of primordial k_i_ would have to be raised to increase 
the pj_ values of the Z'^-boson to the level of the experimental data. Another reason could 
be that soft interactions that fall below the cut-off in p± can build up a p± value of the 
partons before the parton shower starts. 

One specific possibility for too low a rate of emissions in Pythia, explored in this 
thesis, is the following. Pythia version 6.220 uses CTEQ5L parton densities [12]. These 
are leading log, determined with a Aqcd value of 0.192 Gev. However, the parton showers 
in Pythia are not pure leading log. They include corrections to simple DGLAP evolution 
in Q^, that decrease the rate of the x evolution, so that the x evolution is slower than in 
the measurement of the parton densities. This could be compensated by raising the value 
of Aqcd- The higher Aqcd value gives more emissions, which means that a lower value 
of primordial k± is needed for Pythia to give a p^-spectrum similar to the experimental 
one. This could solve the problem with the disagreement between the primordial k± 
value estimated by Heisenberg's uncertainty relation and the experimental value. Which 
corrections Pythia includes is described in Section 4. 

It should also be noted that the primordial kj_ of the initial partons is split between 
the partons in each branching of the parton shower. This means that the partons of the 
hard collision only carry a fraction of the primordial k± of the initial partons 

± J at hard int ~ 1 T I, "-X J initial parton 

(13) 

X'f-initial parton/ 

4 Corrections to pure DGLAP evolution in Pythia 

There are a number of corrections to pure leading log DGLAP backwards evolution in 
Pythia. These, and their effect on the rate of emissions, will be described in this section. 
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Figure 3: Emissions in Pythia are angular ordered: 9i < 62 < 9^. 



4.1 Angular ordering 

Due to coherence effects the opening angles of emitted partons increase with increasing 
virtuality Q^, see Fig. 3. The approximate opening angles can be calculated using x, z, 
and s 

« ^-f^ (14) 

422Q2 + (1 _ z)xH ^ ' 

If an emission has a larger opening angle than the previous one in the Pythia backwards 
evolution, the emission is rejected. This reduces the rate of emissions. 



4.2 The condition that w < 

When one of the quarks in the process qq — > TP emits a gluon this can be seen both 
in a parton shower model, Fig. 4a, and as a 2 — » 2 process. Fig. 4b. The Mandelstam 
variables, Lorentz invariants of the kinematics, can be expressed in terms of the variables 
of the parton shower 



s = (Pl + P2f = ^ (15) 
i = {P,-P,f = -Q^ (16) 

u = {p,-P,)^ = J2m'-s-i = ml-^ + Q^ = Q^-ml^ (17) 

In the picture of a 2 — > 2 process m < 0. For that reason in Pythia there is a condition 
on the emissions that u < 0, using the final expression for u in eq. (17). However, if the 
quark emitting the gluon already before th gluon emission has a negative because 
of preceding shower evolution, the kinematical constraints on u are relaxed and u < 
is just an assumption. No such kinematical constraints are present in the leading log 
DGLAP formalism. This means that the rate of emissions in Pythia is lower than in 
pure DGLAP evolution. 
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Figure 4: The emission of a gluon by one of the quarks at Z°-production seen a) in a 
parton shower model and b) as a 2 ^ 2 process. 




Figure 5: Emission of a gluon by a massless quark. 



4.3 Evolution with as{p\) instead of as{Q^) 

The evolution in Pythia is done with as{p\) instead of as(Q^)- To examine the effect of 
this, an approximate expression for can be derived using the lightcone variables 

p+ = E + p, (18) 



p. = E-p, (19) 

fulfilling 

p+ ■ p- = E"^ - pI = m'^ + p\ (20) 

The kinematics used in Pythia is not this simple but the qualitative result can still be 
derived using lightcone variables. In Fig. 5 these variables are given for the case when a 
massless quark emits a gluon. Conservation of p_ gives 



zp\ (1 — z)p\ zp\_ z{l — z)p^j^ 
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Figure 6: An emitted gluon that develops a final-state shower corresponds to the gluon 
having > 0. 



pl = Q'il-z) (21) 

So < Q^, which means as(pi) > «s(Q^) and thereby more emissions. This is not the 
whole picture, however. As seen in eq. (5) as((5^) has a singularity when = Aqcd- To 
avoid this problem, in Pythia there is a lower cut-off Qq = 1 GeV. The region below 
some cut-off of this order can not be described by perturbative methods. In Pythia this 
region is instead accounted for by non-perturbative methods, such as primordial k±. In 
the case of as{p±) this lower cut-off must be applied to p± instead of Q, so that > Q^. 
As seen in eq. {21) p\ can be below the cut-off for high values of if z is sufficiently close 
to 1. This means that emissions with high values of z are discarded. So the evolution 
with as{p\), instead of as(Q^), has two opposite effects on the rate of emissions. When 
studied the net effect is a decrease, see Subsection 5.1.2. 



4.4 Final-state showers of emitted particles with timelike virtu- 
ality 

Emitted particles that develop final-state showers of their own corresponds to particles 
with timelike virtuality (m^ > for gluons, > for quarks), see Fig. 6 . This does 
not affect the rate of emission but can affect the p_L-spectrum for the produced particle. 
To see this it is again convenient to use lightcone variables. If the gluon in Fig. 5 has a 
positive m^, eq. (21) becomes 

^ _ -Q^+pi m^+p^ _ -Q^ p\ 



zp\ (1 — z)p\ zp\ (1 — z)p\ z{l — z)p\ 

=^ p\ = Q^{\- z)- zm^ (22) 
The introduction of a positive rr? of the gluon clearly causes a decrease in p^. 
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4.5 Other corrections 



There are some other corrections that do not measurably affect either the rate of emissions 
or the p^-spectrum for the produced particle very much. 

The probability of a quark emitting a gluon goes to infinity when the fraction of the 
momentum that the gluon gets in the emission, 1 — z, goes to zero. To avoid this problem, 
in Pythia there is an upper cut-off in z. In the case of a gluon emitting a gluon, the 
singularity occurs both at 2; = and z = 1, and both a lower and an upper cut-off is 
needed. The very soft gluons that fall below this cut-off are resummed. 

Sometimes, e.g. when heavy quarks are produced, the reconstruction of the kinemat- 
ics of a branching fails. The shower is then regenerated. This could affect the rate of 
emissions. When studied however, it does not seem to have any visible effect. 

At large the evolution is matched to first-order matrix elements. This is a separate 
issue from the parton density evolution and will not affect the discussion in this thesis. 

Since quarks carry electric charge they can also emit photons. This is an electromag- 
netic interaction and thereby much less probable than the emission of a gluon. When 
studied, photon emissions do not seem to have any effect on the rate of the x evolution 
either. 

5 Study of Aqcd values, p^-spectra and primordial k± 

First the increase in Aqcd needed to compensate for the corrections to pure DGLAP 
evolution will be quantified in this section. Then the effect of the increase in Aqcd on 
p_L-spcctra will be studied. The hope of this study was to be able to decrease the value of 
primordial k±. This does not seem to be possible, so another approach is taken, the effect 
of final-state showers off emitted particles on p^-spectra. Also a new shower algorithm is 
introduced. 

5.1 Study of the increase in Aqcd needed to compensate for the 
corrections to DGLAP 

5.1.1 Toy model studies 

Some of the main features can be studied in a toy model framework, using only the 
DGLAP equations to generate parton showers. In this framework a certain feature can be 
studied more directly, as opposed to the more comphcated framework of Pythia where 
it is harder to separate the many different effects. This initial study is useful to get some 
understanding of what to expect within a full Pythia simulation. 

The simulated DGLAP evolution should occur at the same pace as the experimentally 
measured parton density scaling violations. Parton density parametrizations are fits to 
the X evolution as a function of Q^, so if the shower does not have the same rate of x 
evolution the right parton density functions would not be obtained. A test of this is to 
pick partons according to the CTEQ5L parton density function at a low Qq, evolve to a 
higher using the toy model and then compare to an x-distribution generated according 
to the CTEQ5L parton density function at the higher Q^. These ^-distributions look close 
to identical, which is comforting. 

The next step is to do backwards evolution. To avoid the problem of quarks originating 
from gluons in the toy model, the backwards evolution is done only for valence quarks. 
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Figure 7: The two x-spectra generated by the toy model, with angular ordering using 
Aqcd = 0.221 GeV and without angular ordering using Aqcd = 0.192 GeV respectively, 
look close to identical. The x-spectra are generated by backwards evolution from the scale 
of the hard coUision, Q"^ — m|, with a center of mass energy of 1800 Gev. For comparison, 
the x-spectrum for the partons at the hard collision is shown. 



where the same quark can be traced through the branchings back to the lower scale Qq. 
The same test as for forward evolution can be done, this time comparing only the x- 
distributions for valence quarks and doing the evolution from the higher scale to the 
lower Qq. Again the result is comforting. 

At this point the effect of angular ordering can be studied. The approximate opening 
angle of each emission is calculated using eq. (14) and an requirement, that the opening 
angle of an emission must be smaller than the previous one in the backward evolution, 
is introduced. This slows down the rate of emissions and consequently the rate of the 
evolution in x. To compensate for this Aqcd can be raised. To quantify the increase in 
Aqcd needed to compensate for angular ordering, x distributions is generated by back- 
wards evolution from an higher scale = m| (to be able to compare to Z° production 
in Pythia) down to a lower scale Qq. The backwards evolution is done with the center 
of mass energy set to 1800 Gev, as is the case for the Tevatron. At first the distribution 
is generated without the angular ordering requirement at the Aqcd value 0.192 Gev used 
in the tuning of the GTEQ5L parton densities, then distributions are generated with the 
angular ordering requirement at different values of Aqcd- The mean values of x can be 
plotted as a function of Aqcd and at some value of Aqcd the mean value of x is the 
same as without angular ordering using Aqcd = 0.192 Gev. According to the toy model 
algorithm, this new Aqcd needed to compensate for angular ordering is 0.221 Gev. As 
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Figure 8: A way of quantifying the increase in Aqcd needed to compensate for angular 
ordering for the process qq Z° at the Tevatron. First the x-spectrum for the initial 
partons is generated without angular ordering using Aqcd = 0.192 GeV. Then x-spectra 
for the initial partons are generated with angular ordering using different values of Aqcd- 
At Aqcd = 0.224 GeV the mean value of x is the same as without angular ordering at 
Aqcd = 0.192 GeV. 



this value is obtained looking only at the mean value of x, it is necessary to check that 
the entire x-spectra look alike. This is done in Fig. 7. 

It should be noted that this study of the increase in Aqcd needed to compensate 
for angular ordering in the toy model only includes valence quarks in the backwards 
evolution. After this initial studies in the toy model, more thorough studies of the effects 
of the corrections to DGLAP evolution on the x evolution has been made using Pythia 
version 6.220, but based on the same principles. 

5.1.2 Pythia studies 

The same way of quantifying the increase in Aqcd needed to compensate for angular 
ordering can be used, with Pythia doing the backwards evolution. If this is done for Z° 
production at the Tevatron the plot in Fig. 8 is obtained. According to Pythia the Aqcd 
needed is 0.224 Gev, not very far from the value given by the toy model. As in Fig. 7 the 
shape of parton densities also agree well. 

This method can be used to find the Aqcd values needed to compensate for all three 
corrections for two different processes, both at the Tevatron and LHC. The Tevatron is a 
pp coUider with a center of mass energy of 1800 Gev and the LHC is a pp collider under 
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construction with a planned center of mass energy of 14000 Gev. It is interesting to see 
how dependent the increase in Aqcd is on the center of mass energy and on the generated 
process. The processes, qq and gg — H°, are chosen to study the differences between 
having qq and gg in the initial state. H° denotes the Higgs boson, which is required within 
the Standard Model but has not yet been found. The Higgs mass is set to 120 Gev, which 
is somewhere in the neighbourhood of where it is expected to be. The results are presented 
in Table 1. 





Tevatron 


LHC 




qq —>■ TP 


gg^RO 


qq Z° 


gg^RO 


angular ordering 


0.224 GeV 


0.239 GeV 


0.221 GeV 


0.236 GeV 


the condition that u < 


0.26 GcV 


0.38 GcV 


0.21 GeV 


0.31 GeV 


as(pi) instead of q;s((5^) 


0.28 GeV 


0.30 GeV 


0.26 GeV 


0.27 GcV 


sum of all effects 


0.31 GeV 


0.50 GeV 


0.29 GeV 


0.46 GeV 



Table 1: The value of Aqcd needed to compensate for different corrections for qq 
and gg H° at the Tevatron and LHC. 

The increase to compensate for angular ordering appears to be almost independent of 
the center of mass energy and not very dependent of which process is generated. It would 
of course be convenient to be able to use the same Aqcd value independently of the center 
of mass energy and particles in the initial state of the process studied. However for the 
other corrections the values seem to fluctuate more. 

Aqcd needs to be increased more when the initial state is gg, compared to having 
the initial state qq. This is not totally unexpected, since quark branchings and gluon 
branchings behave differently. The quark in the branching q — qg tends to take a large 
momentum fraction. This, together with the ordering in Q^, means that the constraints 
on the DGLAP evolution are often fulfilled. In the branching g — > gg one of the gluons 
tends to take a large momentum fraction and the other a small momentum fraction. If 
the gluon of the continued evolution takes a small momentum fraction, the constraints 
are not as naturally fulfilled as for the case with quarks. It would therefore be possible 
to associate different effective Aqcd with different branching processes, but we did not 
pursue it here. There is also a difference in the increase in Aqcd needed between the 
different values of center of mass energy. This difference is not so easy to explain, but 
considerably smaller. 

5.2 The effect of an increased Aqcd value on pj^-spectra and pri- 
mordial /c_L 

In Fig. 9 p_Lz-spectra for the process qq TP at the Tevatron are shown. The p_Lz-spectra 
generated by Pythia, with Aqcd = 0.192 Gev and Aqcd = 0.31 Gev, are compared with 
experimental data [5] . The increase in Aqcd does not seem to affect the p_Lz-spectra much. 
To understand this it could be interesting to look at the effect of Aqcd on at a typical 
scale of the parton shower, = 100 GeV^. When Aqcd is increased from 0.192 Gev 
to 0.31 Gev, aJyQ^ — 100 GeV^) is increased by approximately 10%. The peak position 
of the pxz-spectra is shifted towards higher values by only approximately 1%. This is 
because the p^-kicks caused by the emissions are in different directions in the a;^/-plane 
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Figure 9: p^z-spectra for the process qq TP at the Tevatron, generated by Pythia 
and experimental data. The Pythia predictions are normahzed to the total experimental 
cross section. 



(transverse to the direction of the beam), which means that they partly can cancel each 
other. If you instead look at the Ei_ spectra, the sum of the absolute value of the p_|_-kicks 
caused by emissions, the peak position is indeed shifted towards higher values by the 
expected 10%. 

Formally the p^-kicks corresponds to emission of jets, so in theory the energy of these 
jets could be measured to get an experimental value for E±. Unfortunately, this alternative 
measure of activity in TP events is less accessible experimentally, since also other physics 
could contribute to activity in an event. 

The absent effect of the increased Aqcd value on the p^z-spectra means that something 
else is needed to solve the problem with the too high value of the primordial k±_ needed 
in Pythia to fit experimental data. In Fig. 9 it is also obvioTis that in the p_Lz-spectra 
generated the event rates at the peaks are too high. In Subsection 5.3 another approach 
is described. 

For gg H° at LHC the change in Aqcd has a bigger effect on the pj_H-spectra, see 
Fig. 10. However these spectra is not affected by primordial very much, which can be 
understood as follows. The center of mass energy at LHC is much greater, which means 
more emissions. At each branching the primordial is split between the two particles 
and not much is left at the production of the Higgs particle, cf. eq. (13). 
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Figure 10: p^n-spectra for the process gg — > H° at LHC generated by Pythia. The Higgs 
mass is assumed to be 120 Gev. 



5.3 The effect of final-state showers of emitted particles with 
timelike virtuality on p^-spectra 

As seen in Subsection 4.4 final-state showers of emitted particles with timelike virtuahty 
affect the p^-spectrum of the produced particle. These timelike showers are supposed to 
have virtuality below the spacelike at the emission. It could be argued that this means 
they occur too late to influence the spacelike parton of the branching. The recoil would 
instead be taken by colour- connected neighbour partons. Since both Z° and H*^ are colour 
neutral, their p^-spectra would then be unaffected. If this is right, the p^-spectra can be 
simulated by switching off final-state showers of emitted particles with timelike virtuality. 
A study of the recoiling jets would obviously require more, but we leave that aside here, 
since we do not need to study these jets for the purpose of this thesis. In Fig. 11 this has 
been done and the agreement with experimental data is much better. 

For gg —>■ H*^ at LHC the change in Aqcd has an even bigger effect on the p^-spectra, 
when final-state showers of emitted particles with timelike virtuality are switched off. The 
peak position is then shifted towards higher values by approximately 50%. 

In a new shower algorithm [13] the p^-spectrum of the produced particle would not be 
affected by these timelike showers of emitted particles. This new algorithm is described 
in Subsection 5.4. 
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Figure 11: p_Lz-spectra for the process qq Z° at the Tevatron. Pythia generations with 
and without final-state showers of emitted particles with timelike virtuality, compared to 
experimental data. Both Pythia generations use Aqcd = 0.192 Gev and k± = 2 Gev. 



5.4 p_L-spectra generated by the new algorithm 

In the new algorithm the evolution is done in p']_ = {1 — z)Q'^ instead of Q^. This is 
possible since 

dpi ^ (1 - z)dQ' ^ dQ^ 

pI (1-z)Q2 Q2 i ) 

The ordering of emissions in gives a set of colour connected partons. These partons can 
then evolve a final state shower as a system. The Z°-boson, as a colour neutral particle, 
is unaffected of colour dipoles. 

To compare evolution in p\ with evolution in Q^, the p^z-spectrum generated by the 
new algorithm can be compared to the p^z-spectrum generated by Pythia with final- 
state showers switched off. This is done in Fig. 12. The two spectra look similar except 
for small p±z where the new algorithm has better agreement with experimental data. 

To find which values of Aqcd and primordial k± that give the best fit to experimental 
data using the new algorithm, a test can be performed, is calculated by 

2 ^ " (MC. - data,)^ 
^ ^a2(MQ) + a2(data,) ^ ^ 

where the summation is over the n different experimental data points. In Fig. 13 
values, generated with different values of Aqcd and primordial k±, are shown. The fit 
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Figure 12: The pi^z-spectrum generated by the new algorithm compared to the p±z- 
spectrum generated by Pythia with final-state showers switched ofi^. Both generations 
use Aqcd = 0.192 Gev and k± = 0.4 Gev. 



is best for low values of x^. A perfect fit would give a value equal to the number of 
degrees of freedom (the number of data points minus the number of parameters that are 
to be determined), which in this case is 48. Due to fairly large statistical fiuctuations and 
a lacking description of bin-to-bin correlations some of the values lie below this value. 
The values indicate that again the fit is best for Aqcd ~ 0.19 GeV and k± a; 2.0 GeV, 
but a trend can be seen that increasing Aqcd gives better fits for lower values of k±. 

In Fig. 14 p^z-spectra are generated with some Aqcd and primordial k± values that 
give good values. These p_Lz-spectra are generated using 10 000 000 events, instead of 
the 100000 events used to generate the values in Fig. 13. The values then become 
46.5 for k±^2 Gev, Aqcd = 0.19 Gev, 51.4 for k± = 1.4 Gev, Aqcd = 0.2 Gev and 74.6 for 
k± — 0.6 Gev, Aqcd = 0.22 Gev. In Fig. 14 it can also be seen that the two pj_z-spectra 
generated with k± — 2 Gev, Aqcd = 0.19 Gev and k± — 1.4 Gev, Aqcd = 0.2 Gev 
give good fits to the experimental data, while the one generated with k± = 0.6 Gev, 
Aqcd = 0.22 Gev have too many events at low p±z and too few around the peak. 

Fig. 15 shows the same p_Lz-spectra as in Fig. 14, divided by the experimental data 
points, over a wider range in p±^z- The p_Lz-spectra seem to differ mainly for small values 
of p±z, which is what is expected when different values of primordial k± are used. In 
Fig. 15, there is a trend that the pi^z-spectra with lower values of primordial k± have 
more events with low p±z- A small effect of the change in Aqcd can also be seen for high 
values of p±z in Fig. 15, where higher values of Aqcd mean fewer events around the peak 
and more a bit above that. 
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Figure 13: for fits to experimental data using the new algorithm with different values 
of Aqcd and primordial k± 

At p±z values a bit above the peak many of the generated points lie above the experi- 
mental data points. The fact that the points of the three different generations appear to 
fluctuate in unison owes to the point-to-point fluctuations in the data, where the statistics 
is much smaller than in the Monte Carlo studies. Again the two p^z-spectra generated 
with k± = 2 Gev, Aqcd = 0.19 Gev and k± = 1.4 Gev, Aqcd = 0.2 Gev flt the experi- 
mental data equally well, while the one with k± = 0.6 Gev, Aqcd = 0.22 Gev have too 
many events with low p±z- However it should be noted that the points where the flt of 
the p^z-spectrum with k± = 0.6 Gev, Aqcd = 0.22 Gev is considerably worse than the 
others are mainly the three flrst. A systematic error in the experimental data could affect 
closely lying points in the same way and thus bias the conclusions. 

The new algorithm does not yet include heavy quarks (b, c). When p_L-spectra gener- 
ated by the ordinary Pythia algorithm are studied, the exclusion of heavy quarks does 
not have any measurable effect, however. 

6 Inclusion of heavy quarks in the new algorithm 

This section is an attempt to flnd a way of introducing heavy quarks in the new algorithm. 
No obvious way to do this is found and the work is under development. This is only a 
beginning of these studies. 

6.1 q ^ qg 

To derive the probability for a quark to emit a gluon, the cross section for a process 
that involves a quark emitting a gluon can be divided by the cross section for the same 
process without the emission of the gluon. The process used in this derivation, shown in 
Fig. 16, involves a made-up Higgs particle that couples to the particles needed to make 
the processes as simple as possible with a vertex factor of 1. This is allowed since the aim 
is to derive the probability of the quark emitting a gluon in the colinear limit where it 
does not depend on the rest of the process. 

The needed matrix elements squared are calculated using CompHep [14]. One prob- 
lem of technical character, is that CompHep require the Higgs to be unstable, so we 
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Figure 14: p^z-spectra generated by the new algorithm with Aqcd and primordial kj 
values that give good values. 



have to generate an isotropic Higgs two-body decay. This makes the calculations more 
complicated. 

6.1.1 Derivation of the DGLAP splitting function for massless queirks 

The cross section for the process in Fig. 16b is given by 

a = J dxidx2f{xi)f{x2) as dPSs (25) 

where Xi, X2 arc the fractions of the total proton momentum carried by the initial quarks, 
/ is the structure function of the proton (neglecting the dependence), dPSs is the 
integration over the phase space for the three final particles and the constituent cross 
section is given by 
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^3 oc — — — ()(mi - mn) (26) 
2s 



where |M2_»3p is the matrix element squared and mi is the mass of the intermediate 
particle. 

With the particles denoted according to Fig. 17, the phase space integral is given by 
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Figure 15: The same p_Lz-spectra as in Fig. 14. divided by the experimental data points 
and over a wider range in p±z- The experimental errorbars are also divided by the 
experimental data points. 



where P denotes four-momentum. This can be rewritten as 



1 r c\f 



4 • (27r)2 J ' X1X2S 

where dPS2 is the phase space integral for the process in Fig. 16a and i is the Mandelstam 
variable for the process qq — > H°g. The cross section given in eq. (25) then becomes 



a oc -—— I dxidx2/(xi)/(x2)^^ dmf dPS2 5(mf - ml) (29) 

In the parton shower approximation, the cross section for the process in Fig. 16b can 
also be calculated using the cross section for the process in Fig. 16a and the DGLAP 
splitting function for q — > qg 



aoc / dx,dx,/(xO/(x2)^^d.^gp,_.,,(.)dPS,5(x,x,..-m|) (30) 
In the coUinear limit where 

m? = X1X2ZS =^ dm? = XiX2sdLZ (31) 

eq. (29) becomes 
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Figure 16: a) and b) qq — H° — e+e with and without the emission of a gluon. H° is 
a made-up scalar particle that has the couplings necessary for the processes. 




Figure 17: The final particles of the process in Fig. 16b. 3 is the gluon, i is the intermediate 
made up Higgs boson, 1 and 2 are the e"'"e~-pair and denotes the available energy or 
momenta from the collision. 
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(32) 



4- (2^)2 

The condition that the two expressions for the cross section, eq. (30) and eq. (32), should 
be equal gives, after some further considering 



P. 



q^qg 



Q' |M2_ 



(33) 



2^2 



The matrix elements squared can be calculated, e.g. using CompHep. The process in 
Fig. 16b can not be separated from the process where the other quark emits the gluon. 
The matrix elements must be added and then squared. In the collinear limit where the 
gluon is emitted in the direction of the first quark the matrix element squared of the 
second process disappears. The interference term however can contribute to the derived 
expression for the splitting function. Using the CompHep matrix elements squared gives 
the following expression for the splitting function 



q^qg 



4 s^{z^ + 1) - 



(34) 



3s(s(l-z) + Q2) 

In the collinear limit, where is small compared to s, eq. (34) goes to the DGLAP 
splitting function in eq. (2). 
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6.1.2 Massive quarks 



The corresponding calculation can be done for the case when the quark emitting the gluon 
in the process in Fig. 16b has a non- negligible mass. For a c quark the calculation gives 

n^cgU) - 3 (^2 + g2)2((l_^)^_^2_Q2) 

^ dg2 ~ 27r3 (m2 + Q2)2((i_^)^_^2_g2) I ) 

Since the new algorithm is ordered in p^, what is needed is really ^'^pT"'^ instead of '^^Iq-i^ ■ 
For the case of a c quark emitting a gluon as a function of becomes 



((s-ml)(l-z)-Q')Q^s 
{§ - mlY 

g2 ^ V |V ; _ _ ,^2)^ _(1 -zY-i^ 

dQ2 1 ,^ _ „;2 



dpi 2 
Neglecting the term gives 



1(1--)^-^ 



(36) 



(s — m^)pj 



{l-z)s 
dQ^ s — m 



dp2 (1 — z)s 



(37) 



Even with these simplified expressions for Q"^ and the expression for ^^"^ gets very 
complicated. 

Complexity in itself is not an absolute hindrance, however: if a simple upper estimate 
can be found, standard rejection techniques can be used to obtain results in agreement 
with the full expression. It therefore remains to find some suitable such estimate in a 
p_L-oriented framework. 



6.2 g ^ qq 

To derive the probability for a gluon to split into a cc pair, the same method as in 
Subsection 6.1 can be used, this time with the processes shown in Fig. 18. The following 
expression is obtained when, again using CompHep, the matrix element squared for the 
process in Fig. 18b is divided by the matrix element squared for the process in Fig. 18a 



4770;.= 



{ml + Q2)2 _ ^2) 



+ 47rQ; 



mi 



2zQ'^ - 2m2Q2 _ _ 2m\ - 2m\z 



{ml + Q2)2 _ ^2) 



+ 



(38) 
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Figure 18: qq — H^ ^ e^e . In b) the quark comes from a gluon splitting into a qq-pair. 
H" is again a made-up scalar particle that has the coupling necessary for the processes. 

If the fact that there is a massive c quark in the initial state of the process in Fig. 18a 
is taken into account the denominator in Eq. 26 should be 2{s — ml) instead of 2s. This 
means that the expression in Eq. 38 should be multiplied with ^ ^j^ich gives the 

following expression 



Ana. 



Q\z' + {l-z) 



+ (l + 2z^W 



(ml + Q2)2 
1 - z)Q^ + ml {ml + Q^) + m^z 



+ 



(39) 



where the first term takes the form of the DGLAP splitting function when ml is set to 0. 

In the process in Fig. 18b there is a massive c quark in the final state. This means 
that the phase space integral gets more complicated 



dPS, = dPS, 



4- (2^)2 



{xiX2S — mf + m2)2 — AxiX2sml 
X1X2S — mf + ml 



dm. 



dt 



X1X2S 



(40) 



Again it remains to find a simple upper estimate as a starting point for a more precise 
inclusion of mass effects. 



7 Summary 

To explore perturbative and nonperturbative effects in transverse momentum generation 
P±-spectra have been studied. More specificaly, the influence of corrections in Pythia to 
pure leading log DGLAP evolution on x evolution and generation of p^-spectra has been 
investigated. This has been done in hope that the reason for the need of a primordial k± 
much larger than expected might be found. 

To compensate for the corrections to pure leading log DGLAP evolution in Pythia, 
Aqcd needs to be raised to the values given in Table 1. For qq at the Tevatron 

this increase in Aqcd does not affect the generated p_Lz-spectra very much. This can be 
understood as the effect of kicks of the branchings in different directions in the plane 
transverse to the beam direction partly canceling each other. The generated p^z-spectra 
do not fit the experimental data very well. For gg H^ at LHC a changed Aqcd has 
a larger effect on the generated p^n-spectra, due to higher center of mass energy and 
consequently more branchings. In this case the primordial kx. of the initial partons does 
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not affect tfie generated p^n-spectra very mucli, since it is split between tfie partons of 
the many branchings. 

Final-state showers of emitted particles with timelike virtuality affect the generated 
Px-spectra. It could be argued that these showers occur too late to influence the spacelike 
parton of the branching, since they are supposed to have virtuality below the virtuality 
of the branching. When these are turned off the fit to the experimental data for qq — > TP 
at the Tevatron is much better. 

In a new shower algorithm the branchings are ordered in instead of . This gives 
a set of colour connected partons which can then undergo final state showers as a system. 
These final state showers would not affect the produced Z", since it is a colour neutral 
particle. In this algorithm the increase in Aqcd has a larger effect on the generated 
P±z-spectra. When three different p^z-spectra, using values of Aqcd and primordial 
that give good values, arc studied the fit is equally good for the two spectra using 

= 2 Gcv, Aqcd = 0.19 Gev and = 1.4 Gev, Aqcd = 0.2 Gev respectively, while the 
fit is a bit, but not very much, worse for the one using k± = 0.6 Gev, Aqcd = 0.22 Gev. 
The problem with a large value of primordial k± is not completely solved but reduced. 

The studies to include heavy quarks in the new algorithm are under development. No 
obvious way of doing this has yet been found. 



26 



References 

[1] T. Sjostrand, P. Eden, C. Friberg, L. Lonnblad, G. Miu, S. Mrenna and E. Norrbin, 
Comput. Phys. Commun. 135 (2001) 238 [arXiv:hep-ph/0010017]. 

[2] T. Sjostrand, L. Lonnblad, S. Mrenna and P. Skands, arXiv:hep-ph/0308153. 

[3] F. Abe et al. [CDF Collaboration], Phys. Rev. Lett. 66 (1991) 2951. 

[4] B. Abbott et al. [DO Collaboration], Phys. Rev. Lett. 80 (1998) 5498 [arXiv:hep- 
ex/9803003]. 

[5] T. Affolder et al. [CDF Collaboration], Phys. Rev. Lett. 84 (2000) 845 [arXiv:hep- 
ex/0001021]. 

[6] G. AltarelU and G. Parisi, Nucl. Phys. B126 (1977) 298; 

V.N. Gribov and L.N. Lipatov, Sov. J. Phys. 15 (1972) 438, ibid. 75; 
Yu. L. Dokshitzer, Sov. J. Phys. JETP 46 (1977) 641 

[7] B. R. Webber, Ann. Rev. Nucl. Part. Sci. 36 (1986) 253. 

[8] T. Sjostrand, Phys. Lett. B 157 (1985) 321. 

[9] M. Bengtsson, T. Sjostrand and M. van Zijl, Z. Phys. C 32 (1986) 67. 

[10] G. Miu and T. Sjostrand, Phys. Lett. B 449 (1999) 313 [arXiv:hep-ph/9812455]. 

[11] C. Balazs, J. Huston and I. Puljak, Phys. Rev. D 63 (2001) 014021 [arXivrhep- 
ph/0002032]. 

[12] H. L. Lai et al, Eur. Phys. J. C 12 (2000) 375 

[13] T. Sjostrand, LU TP 04-05 [hep-ph/0401061] and in preparation 

[14] A. Pukhov et al, arXiv:hep-ph/9908288. 



27 



hep-ph/0401121 LU TP 04 - 01 
January , 2004 



Perturbative and Nonperturbative 
Effects in Transverse Momentum 

Generation 

Master of Science Thesis by Erik Thome 

Thesis advisor: Torbjorn Sjostrand 

Department of Theoretical Physics, 
Lund University, Lund, Sweden 



Abstract 

The transverse momentum of a colour-singlet massive particle in a hadronic 
collision is built up by two components, the perturbative effect of parton 
branchings and the nonperturbative effect of primordial k±. In previous 
studies of transverse momentum spectra for Z° production at the Tevatron, 
the best fit to the experimental data are given when the primordial k± is set 
to a much higher value than what is expected considering the confinement 
of partons in the proton. We here investigate the possibility that the reason 
for this is that too few branchings are generated in showers, compared to 
the evolution used in the tuning of parton densities. This could then be 
compensated by increasing the value of Aqcd- The study is done using 
the regular Pythia showering routines and a new algorithm where the 
branchings are ordered in transverse momentum instead of virtuality Q^. 



1 Introduction 



Again, of bodies some are composite, others the elements of which these 
composite bodies are made up. These elements are indivisible and unchange- 
able, and necessarily so, if things are not all to be destroyed and pass into 
non-existence, but are to be strong enough to endure when the composite bod- 
ies are broken up, because they possess a solid nature and are incapable of 
being anywhere or anyhow dissolved. It follows that the first beginnings must 
be indivisible, corporeal entities. 

Epicurus "Letter to Herodotus" 
approximately 300 B.C. 

What are the building blocks of the universe? That is a very profound question, 
pondered upon by philosophers through the milleniums. For a long time the theories were 
purely metaphysical, with no possibility whatsoever to be tested experimentally. Thanks 
to major breakthroughs in theoretical and experimental physics during the twentieth 
century we now have the Standard Model, describing those elementary particles and their 
interactions. This model has been extremely successful in predicting the outcome of 
experiments, performed at increasingly large accelerators. 

The elementary particles are divided into quarks, leptons and gauge bosons. Also the 
corresponding antiparticles exist. The gauge bosons mediate the four fundamental forces 
of nature: electromagnetism, the weak force, the strong force and gravity. The interesting 
particles for this thesis are the quarks and the gluons, the gauge bosons mediating the 
strong force. These are denoted q and g, with q denoting an antiquark. The theory 
describing the interactions of quarks and gluons is called quantum chromo dynamics, 
abbreviated QCD. 

Free quarks have never been observed. The theoretical explanation for this has to 
do with the behaviour of the strong force. Here, it is useful to make a comparison with 
electromagnetism. Instead of electromagnetic charge, the strong force acts on particles 
with colour charge. This is, of course, not real colour, but is named this way since the 
behaviour of the colour charges resembles the behaviour of real colours in the way that 
mixing all the colour charges makes a colour neutral object. Of the elementary particles 
it is the quarks and gluons that carry colour charge. The crucial point is that the force 
between two colour charged particles does not decrease with the distance separating them, 
as is the case for electromagnetism. This means that it takes an infinite amount of energy 
to separate a qq pair, if they were to be connected to each other by a colour field. When 
the distance between the quarks gets sufficiently large, the potential energy becomes large 
enough to create a new qq pair. The new qq pair screens the colour field between the 
original charges. But now there are two colour connected qq pairs. So when the quark 
and the antiquark of a colour connected qq pair are moved away from each other, new qq 
pairs are created, but free quarks will not be obtained. 

Another important difference from electromagnetism is that the particle mediating the 
strong force, the gluon, is not a colour neutral object. This means that a gluon can emit 
a gluon. So there are three different QCD branchings of particles: q — > qg and g — > qq, 
which have their corresponding branchings in electromagnetism (q — > q7 and 7 — > qq 
respectively), and g — > gg, which has no similar branching in electromagnetism. 

Quarks can carry the colour charges red, green or blue, the antiquarks the complemen- 
tary anticolours and gluons one colour and one anticolour. As noted above only colour 
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neutral objects can appear as free particles. Colour neutral particles out of quarks can be 
obtained in two different ways, either qq (colour-anticolour) or qqq (the mixture of red, 
green and blue makes a colour neutral object). The former particles are called mesons 
and the latter baryons (of course also antibaryons qqq exist), of which the proton and 
the neutron are the most well-known. The common name for all particles consisting of 
quarks is hadrons. 

There are three generations of quarks ^ j ^ 5 ( ^ ) ^ ^ ^ ' ^' ^ ^ 

quarks are fairly light and will be assumed to be massless in this thesis, the c and b 
quarks are considerably heavier and the t quark is very heavy and too shortlived to be 
bound in any hadron. 

To test the Standard Model experimentally, particles are accelerated to high energies 
and then collided. In the collisions the energy can be used to create new particles. In this 
way many particles predicted by the Standard Model have been found. Much of the effort 
in particle physics today is aimed at finding the Higgs particle, the last particle predicted 
by the Standard Model that has not yet been found. This search has been performed 
at the Tevatron, an accelerator at Fermilab outside of Chicago, and will be continued at 
LHC (Large Hadron Collider), an accelerator under construction at CERN in Geneva. 

At the Tevatron protons and antiprotons are collided and at LHC protons are collided 
with protons. The energies are high enough to resolve the quarks and the gluons in the 
proton; these are called partons with a common name. The process of interest is often 
a hard collision between partons, which is not easy to extract from the experimental 
data or calculate theoretically. Firstly, in the final state, the partons created at the 
high energy scale of the hard collision branch into showers of partons. These parton 
showers evolve to a lower energy scale, where the partons are bound into hadrons. It is 
these hadrons that show up in the detector. The hadronization of partons can not be 
described by perturbative methods. Secondly, in the initial state, when the accelerated 
particles are protons, there is a problem of knowing the inner structure of the proton at 
different energy scales. High energy scales correspond to short distance interactions, as 
hard collisions of partons. Since the process of interest is a hard collision of partons, the 
distribution of partons in the proton at this energy scale must be known. The partons of 
the hard collision can be generated by letting the partons of the proton at a lower energy 
scale develop initial-state parton showers up to the scale of the hard collision. These 
initial-state parton showers will be the subject of the study in this thesis. 

At the scale of particle physics our universe is described by quantum mechanics, which 
means that it is not possible to predict the outcome of a single event. Only the probabili- 
ties of what is going to happen can be predicted. It is therefore necessary to average over 
a large number of events in a meaningful way. Perturbative methods can be applied to 
describe events down to a low energy scale, under which nonperturbative methods have 
to be used. In an event hundreds of particles can be produced. This means that even if 
the calculations could be done using perturbative methods, the complexity would exclude 
many approaches. One approach that does work is doing Monte Carlo simulations of a 
large number of events. Both initial-state and final-state parton showers are generated. 
This way spectra of different variables, that can be compared to experiment, are obtained. 
Fvent generators are used to interpret the outcome of experiments and to predict what 
to look for in the future. One of the most successful event generators, that will be used 
in this thesis, is called Pythia [1, 2]. 

One interesting variable to simulate in event generators is transverse momenta p_L, the 
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momentum that particles have transverse to the accelerated beam. In the branchings the 
partons undergo, they get recoil transverse to the beam axis. This means that the more 
branchings, the more p± of the partons. However there is also an opposite effect, the 
accumulated p± of a parton is split between the two new partons of a branching. Still, 
the branchings affect so to look at p± is a good way to find out something about the 
branchings of an event. p± can be seen as a measure of how violent a process is. 

One process where it is possible to make a comparison between theory, event generators 
and experimental data in a rather clean way is qq — > Z°. When the Z° decays to an e"'"e~ 
pair, its p± is reasonably easy to measure. The p± of the Z° boson is the sum of the 
p± of the partons in the hard collision, where it was created. The p± of these partons 
is built up by two components; a perturbative component from the parton branchings 
and a non-perturbative component called primordial k± to be described as follows. The 
partons in the proton have an uncertainty in momenta that, according to Heisenberg's 
uncertainty relation, should be inversely proportional to the radius of the proton. This 
gives rise to a transverse momentum, that is called primordial k±. For high p± values of 
the Z° boson the perturbative component from the parton branchings dominates, while 
the low p± part of the p_L-spectrum is sensitive to primordial k±. 

The pj_ of the Z° and bosons in the process qq — > Z°/W^ have been measured 
at the Tevatron [3, 4, 5] and there is a problem with the Pythia prediction of the p±- 
spectrum for small p± as follows. To get the best fit to experimental data for small p± the 
value of the primordial k± has to be set considerably higher than expected, considering 
the radius of the proton. 

This thesis investigates one possible cause to this problem. The study will be made 
for the case of Z° production. In Pythia, some perturbative branchings, that normally 
are assumed to contribute, are not simulated. This means that the p_L-spectrum of Z° 
boson will be shifted to lower values of p±. To compensate for this the other component 
building up the p± of the Z° boson, the primordial A;_l, has to be increased. 

If the cause of the problem is that Pythia generates too few branchings, there might 
be a way of solving it. In the expression for the running strong coupling constant as{Q^) 
there is a parameter Aqcd that has to be measured experimentally. An increase of Aqcd 
increases the rate of branchings. At first sight this does not seem to provide any solution 
to the problem, since the same Aqcd used to determine the parton densities, which are 
used in Pythia, must be used in the Pythia parton showers. However, the parton 
densities are leading log. In Pythia there are some corrections to pure leading log 
DGLAP evolution, introduced to give better agreement with next to leading log results. 
It is these corrections that decrease the rate of branchings, so to compensate for this Aqcd 
can be increased. What is studied in this thesis is how much Aqcd should be increased to 
compensate for the corrections and to what extent this solves the problem with primordial 
k^. 

In Section 2 the concepts of parton showers and the DGLAP equations are introduced 
and a more specific description of the way of doing initial-state parton showers in Pythia 
is given. Section 3 is a description of the problem with the too high value of primordial 
A;j_, when comparing to experiment. The corrections to DGLAP evolution in Pythia are 
described in Section 4. The study of the generation of p_L-spectra is presented in Section 5. 
In this section also a new algorithm for doing initial-state showers is introduced. This 
new algorithm does not include heavy quarks. Section 6 is a attempt to find a way of 
introducing heavy quarks in the new algorithm. And finally Section 7 is a summary of 
the results of this thesis. 
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Figure 1: From the three processes q — > qg, g — > gg and g — > qq it is possible to build up 
complicated showers. The emissions are ordered in . 

2 Parton showers 

To simulate events, that can be compared with the data from accelerators, either matrix 
element calculations or parton showers can be used. In the matrix element approach 
amplitudes and phases are used to do a fully quantum mechanical calculation. Calcu- 
lating matrix elements order by order in the strong coupling constant 0,^ is a good, but 
complicated, method to get cross sections for processes. As each branching contributes 
with one order in ag, the matrix element calculations get too complicated for processes 
involving more than a few branchings. An approximate approach is needed. In such an 
approach probability distributions, instead of amplitudes and phases, can be used. This 
is an ideal method for doing computer simulations. Using the probabilities for the differ- 
ent branchings, given by the DGLAP equations [6], one parton branching into a shower 
of partons can be simulated [7]. This method can be applied to events with arbitrary 
many branchings, but is only accurate to what is called leading log, where to each order 
in ftg, only the most divergent contribution is described correctly. In this parton shower 
approach the histories of the partons are given, so they can be traced back through the 
branchings. 

2.1 General concepts of parton showers 

In the introduction the proton is said to consist of three quarks; this is a simplification. 
In the proton different processes is constantly occuring: q — > qg, g — > gg, g — > qq and 
their inverses. The quantum numbers of a quark originating from a gluon splitting are 
always compensated by the quantum numbers of the antiquark originating from the same 
splitting. These quarks are called sea quarks, while the three quarks that determine the 
quantum numbers of the proton are called valence quarks in analogy with atomic physics. 

As a consequence of these constantly occuring processes, the structure of the proton 
depends on the momentum scale it is probed by. A higher can resolve more gluons 
and sea quarks. To describe the structure of the proton, parton density functions fa{x^ Q^) 
are used. They express the probability of finding a parton of type a with a fraction x 
of the total proton momentum in a proton probed with Q^. To calculate parton density 
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functions nonperturbative input is needed. This nonperturbative input is equivalent to 
having the fa{x^Q^) defined at one = Qq. Given such initial conditions, the parton 
density function at another can be calculated, see Subsection 2.2. The dependance 
of the parton density functions is often called scaling violations, since the are no longer 
functions of x alone. 

Another way of describing this gradual resolution of the proton is to start with one 
initial parton and do an evolution in Q^. As increases the parton branches into more 
and more partons and a shower of partons arises. Such a parton shower is built up by 
the three processes (if processes involving photons are neglected) q — > qg, g — > gg and 
g — > qq, see Fig. 1. To do a simulation of parton showers, the probability for these 
processes as a function of is needed. In Subsection 2.2 the DGLAP equations, giving 
these probabilities, are introduced. 

This parton shower approach can be used to generate events at particle colliders. An 
event is separated into one hard process, with the highest Q^, and many softer interactions 
both before and after the hard process. To describe these softer interactions, the parton 
shower approach is useful. The parton showers of an event can be separated into final-state 
showers and initial-state showers. The final-state showers evolve if the final particles of the 
collision are quarks or gluons. This can be simulated by an evolution in Q^, starting with 
the quarks at the scale of the hard collision and going down to the scale of hadronization. 
If the colliding particles are protons there is also initial-state parton showers, starting 
with the partons at some low scale Qq going up to the scale of the hard collision. In Fig. 1 
the parton shower starts with a quark 1 at Qo, which branches into a quark 2 and a gluon 
3 at a higher Q^. At an even higher quark 2 branches into a quark 4 and gluon 5 and 
at another gluon 3 splits into a quark 6 and an antiquark 7. The branchings continue 
up to the scale of the hard collision. Only one of the partons in the cascade takes part 
in the hard collision. When the partons that participates in the hard collisions are given, 
they can be traced back to the low scale by backwards evolution, to be described below. 
The backwards evolution starts with the parton of the hard collision, for instance gluon 
8, tracing it to the branching of gluon 5 at a lower Q^. Gluon 5 is emitted from quark 
2 at an even lower Q^. Quark 2 originates from the branching of quark 1, which is the 
initial parton of the shower at scale Q^. In this way the p± of the colliding particles, and 
consequently the p± of the produced particle, can be calculated. 



In DGLAP evolution the emissions are ordered in spacelike virtualities and the prob- 
ability of the process a ^ be is given by 



where Pa^bc{z) is the splitting function depending on the fraction of the momentum of 
particle a taken by particle b. For the processes that are used to build up a QCD shower, 
the splitting functions are given by 



2.2 DGLAP equations 







(3) 



5 



^g^-) = \{z' + {l-zf) (4) 
In eq. (1) the running strong coupling constant is used 

where nf is the number of flavours and Aqcd is a parameter of the Standard Model that 
has to be measured experimentally. Aqcd will be of great importance for this thesis. 

The parton shower approach is an exclusive approach where one parton, that branches 
into a shower as increases, is studied. To include all the partons of the proton, one can 
study the behaviour of parton densities as increases. Beginning with initial conditions 
for the parton densities at one Qo, considering the evolution of all the partons, the parton 
densities at another can be calculated. The dependence on of the parton density 
for a particle of type h at the value x is given by 

2.3 Initial-state showers in Pythia 

In this thesis initial-state showers in Pythia [8, 9, 10] are studied. Pythia uses the 
parton shower approach to do Monte Carlo simulations of events. The DGLAP equations 
are used to simulate at which Q^, and with which the branchings occur. 

In the evolution a branching can only occur if the parton has not already branched. 
This introduces a Sudakov form factor to the expression for the probability 

dn^.. = -qT^^P^^^^(^) [- ^ j (7) 

To do a Monte Carlo simulation, a function N[Q'^) for the number of partons that has 
not branched is needed 



N{Q') = N{Qt>){l - dV) 



This number is then set to a random number R times the number of partons at Qq, 
N[Qq). So the at which the branching occurs is given by 



dV = 1- R = R 



using the fact that 1 — i? is a random number between and 1 if i? is it. The z of the 
branching can then be generated by setting 

r dz'Pa^,,{z') = rC dz'Pa^Uz') (10) 

Jo Jo 
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As seen in eq. (2) and eq. (3) the probability of emissions of very soft gluons, where the 
gluon takes a very small part of the momentum, goes to infinity. It is therefore necessary 
to introduce a cut-off e in z. In the case of q — > qg only an upper cut-off is needed, while 
in the case of g — > gg both an upper and a lower cut-off is needed. Instead of integrating 
from on the left side and up to 1 on the right side in eq. (10), the integration is from e 
and up to 1 - e respectively. 

In the forward evolution is generated this way for the different possible branchings 
(q — > qg, g — > gg and g — > qq) and the one with the lowest is chosen. Forward evolution 
can be very inefficient when a specific kind of hard process is studied, since the final-state 
partons are not known beforehand. If no final-state partons with the right properties 
to participate in the hard process are generated in the evolution, it can not be used. A 
more efficient method, used in Pythia, is to fix the hard scattering using evolved parton 
densities and thereafter use backwards evolution to form this into exclusive events. Such 
an approach can again be formulated in Monte Carlo terms, but now starting from eq. (6) 
instead of eq. (1). This gives 

where z has been rewritten as x, the momentum fraction of the parton before (in the 
ordering of the backwards evolution) the branching, divided by y, the momentum fraction 
of the parton after the branching. In the exponential the probability is integrated from 

to Q^ax' virtuality of the preceding step, instead of from Q\ to Q^, as is the case 
for forward evolution. 

The same procedure can then be repeated starting at the of the branching. In 
this way the partons of the hard interaction are evolved back to the initial partons at the 
lower scale Q\. 



3 Primordial 

The partons are confined in the proton of radius 0.7 fm. According to Heisenberg's 

uncertainty relation this means that they have a momentum inside the proton. The part 
of this momentum that is in the beam direction gives the fraction x of the total proton 
momentum carried by the parton, described by parton density functions. But there is 
also a part of this momentum that is transverse to the beam direction. This transverse 
momentum of the partons, due to the confinement in the proton, is called primordial A;_l 
and its mean value can be estimated by 

, h 0.2GeV-fm ^ , , 

IkA^ — ^ — ^ 0.3 GeV 12 

Tp 0.71m 

To compare this value to experiment, we can look at the p_L-spectrum of the produced 
Z° boson in the process qq — > Z°. The reason for looking at this particular process is 
that when the Z° boson decays to a e~ pair the p_L of the Z° boson is reasonably easy 
to measure. The p_L-spectrum of the Z° boson depends on the value of the primordial 
A;_L and on the p_L-kicks the quarks get when they emit gluons. For high p_L values of 
the Z° boson the effect of the p_L-kicks in the shower dominates, while the low p_L part 
of the p_L-spectrum is sensitive to primordial k^. Fig. 2 gives a schematical picture of 
this. In one set of data from the Tevatron the Pythia prediction of the pj_-spectrum 
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Primordial k 



_L 




Primordial k± \ 9 ' " ^~"^'"'*h"*i»«Jf 

Figure 2: Production of a Z°-boson in a parton shower model. At each branching the 
accumulated p_L is split between the particles in z vs. 1 — z fractions. New p± values 
for the partons after the branching are generated, and the shower continues. This means 
that showers with many branchings are not very sensitive to the value of primordial k±. 

fits the experimental data best for a primordial k± value of 2.15 Gev [11]. This value is 
considerably higher than the value estimated by Heisenberg's uncertainty relation. One 
reason for this could be that Pythia has too low a rate of emissions. Then there would 
be too few p_L-kicks and the value of primordial k± would have to be raised to increase 
the p± values of the Z°-boson to the level of the experimental data. Another reason could 
be that soft interactions that fall below the cut-off in p± can build up a p± value of the 
partons before the parton shower starts. 

One specific possibility for too low a rate of emissions in Pythia, explored in this 
thesis, is the following. Pythia version 6.220 uses CTEQ5L parton densities [12]. These 
are leading log, determined with a Aqcd value of 0.192 Gev. However, the parton showers 
in Pythia are not pure leading log. They include corrections to simple DGLAP evolution 
in Q^, that decrease the rate of the x evolution, so that the x evolution is slower than in 
the measurement of the parton densities. This could be compensated by raising the value 
of Aqcd- The higher Aqcd value gives more emissions, which means that a lower value 
of primordial k± is needed for Pythia to give a pj_-spectrum similar to the experimental 
one. This could solve the problem with the disagreement between the primordial k± 
value estimated by Heisenberg's uncertainty relation and the experimental value. Which 
corrections Pythia includes is described in Section 4. 

It should also be noted that the primordial k± of the initial partons is split between 
the partons in each branching of the parton shower. This means that the partons of the 
hard collision only carry a fraction of the primordial k± of the initial partons 



4 Corrections to pure DGLAP evolution in Pythia 

There are a number of corrections to pure leading log DGLAP backwards evolution in 
Pythia. These, and their effect on the rate of emissions, will be described in this section. 




(13) 
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Figure 3: Emissions in Pythia are angular ordered: 6i < 62 < 03- 



4.1 Angular ordering 

Due to coherence effects the opening angles of emitted partons increase with increasing 
virtuality Q^, see Fig. 3. The approximate opening angles can be calculated using x, z, 
and s 

^ Az-^Q"^ + z)x'^s ^ ' 

If an emission has a larger opening angle than the previous one in the Pythia backwards 
evolution, the emission is rejected. This reduces the rate of emissions. 



4.2 The condition that m < 

When one of the quarks in the process qq — > Z° emits a gluon this can be seen both 
in a parton shower model, Fig. 4a, and as a 2 — > 2 process. Fig. 4b. The Mandelstam 
variables, Lorentz invariants of the kinematics, can be expressed in terms of the variables 
of the parton shower 



s = {P, + p,f = !!^ (15) 

z 

i = {P,-Psf = -Q' (16) 
u = (p,-P^Y = Y^m'-s-i=ml-'^ + Q' = Q'-ml^-^ (17) 

In the picture of a 2 — > 2 process u < 0. For that reason in Pythia there is a condition 
on the emissions that u < 0, using the final expression for ii in eq. (17). However, if the 
quark emitting the gluon already before th gluon emission has a negative because 
of preceding shower evolution, the kinematical constraints on ii are relaxed and u < 
is just an assumption. No such kinematical constraints are present in the leading log 
DGLAP formalism. This means that the rate of emissions in Pythia is lower than in 
pure DGLAP evolution. 
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Figure 4: The emission of a gluon by one of the quarks at Z°-production seen a) in a 
parton shower model and b) as a 2 — > 2 process. 




Figure 5: Fmission of a gluon by a massless quark. 



4.3 Evolution with asip^) instead of as{Q'^) 

The evolution in Pythia is done with as{p']_) instead of as{Q^). To examine the effect of 
this, an approximate expression for p\ can be derived using the lightcone variables 

P+ = E + p, (18) 



P- = E-p, (19) 

fulfilling 

pj, ■ p_ = E"^ - pI = + pI (20) 

The kinematics used in Pythia is not this simple but the qualitative result can still be 
derived using lightcone variables. In Fig. 5 these variables are given for the case when a 
massless quark emits a gluon. Conservation of p_ gives 

, ^ -Q' + pI I pI _ -Q' I pI 



zp^ (1 — z)p^ zp^ zil — z)p\ 
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Figure 6: An emitted gluon that develops a final-state shower corresponds to the gluon 
having > 0. 



pI = Q'{1-z) (21) 

So p\ < Q^^ which means as{p']_) > as{Q^) and thereby more emissions. This is not the 
whole picture, however. As seen in eq. (5) as{Q^) has a singularity when = Aqcd- To 
avoid this problem, in Pythia there is a lower cut-off Qq = 1 GeV. The region below 
some cut-off of this order can not be described by perturbative methods. In Pythia this 
region is instead accounted for by non-perturbative methods, such as primordial k±. In 
the case of o.s[p\) this lower cut-off must be applied to pj_ instead of so that p^j_ > Qq. 
As seen in eq. (21) p\ can be below the cut-off for high values of if z is sufficiently close 
to 1. This means that emissions with high values of z are discarded. So the evolution 
with as{p']_), instead of as{Q^), has two opposite effects on the rate of emissions. When 
studied the net effect is a decrease, see Subsection 5.1.2. 



4.4 Final-state showers of emitted particles with timelike virtu- 
ality 

Emitted particles that develop final-state showers of their own corresponds to particles 
with timelike virtuality (m^ > for gluons, > for quarks), see Fig. 6 . This does 
not affect the rate of emission but can affect the p_L-spectrum for the produced particle. 
To see this it is again convenient to use lightcone variables. If the gluon in Fig. 5 has a 
positive m^, eq. (21) becomes 

^ _ -Q^ +pI + p\ _ -Q^ pI 

"•"/-I \n n~r/-i \n~r 



zpl (1 - z)p1 zp\ (1 - z)p\ z{l - z)p1 



^ pI = Q\1 -z)-zm^ (22) 
The introduction of a positive of the gluon clearly causes a decrease in p^. 
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4.5 Other corrections 



There are some other corrections that do not measurably affect either the rate of emissions 
or the p_L-spectrum for the produced particle very much. 

The probability of a quark emitting a gluon goes to infinity when the fraction of the 
momentum that the gluon gets in the emission, I — goes to zero. To avoid this problem, 
in Pythia there is an upper cut-off in z. In the case of a gluon emitting a gluon, the 
singularity occurs both at 2; = and 2; = 1, and both a lower and an upper cut-off is 
needed. The very soft gluons that fall below this cut-off are resummed. 

Sometimes, e.g. when heavy quarks are produced, the reconstruction of the kinemat- 
ics of a branching fails. The shower is then regenerated. This could affect the rate of 
emissions. When studied however, it does not seem to have any visible effect. 

At large the evolution is matched to first-order matrix elements. This is a separate 
issue from the parton density evolution and will not affect the discussion in this thesis. 

Since quarks carry electric charge they can also emit photons. This is an electromag- 
netic interaction and thereby much less probable than the emission of a gluon. When 
studied, photon emissions do not seem to have any effect on the rate of the x evolution 
either. 

5 Study of Aqcd values, ]9x-spectra and primordial k_i 

First the increase in Aqcd needed to compensate for the corrections to pure DGLAP 
evolution will be quantified in this section. Then the effect of the mcrease m Aqcd oh 
p_L-spectra will be studied. The hope of this study was to be able to decrease the value of 
primordial k±. This does not seem to be possible, so another approach is taken, the effect 
of final-state showers off emitted particles on p_L-spectra. Also a new shower algorithm is 
introduced. 

5.1 Study of the increase in Aqcd needed to compensate for the 
corrections to DGLAP 

5.1.1 Toy model studies 

Some of the main features can be studied in a toy model framework, using only the 
DGLAP equations to generate parton showers. In this framework a certain feature can be 
studied more directly, as opposed to the more complicated framework of Pythia where 
it is harder to separate the many different effects. This initial study is useful to get some 
understanding of what to expect within a full Pythia simulation. 

The simulated DGLAP evolution should occur at the same pace as the experimentally 
measured parton density scaling violations. Parton density parametrizations are fits to 
the X evolution as a function of Q^, so if the shower does not have the same rate of x 
evolution the right parton density functions would not be obtained. A test of this is to 
pick partons according to the CTEQ5L parton density function at a low Qq, evolve to a 
higher using the toy model and then compare to an x-distribution generated according 
to the CTEQ5L parton density function at the higher Q^. These ^-distributions look close 
to identical, which is comforting. 

The next step is to do backwards evolution. To avoid the problem of quarks originating 
from gluons in the toy model, the backwards evolution is done only for valence quarks. 



12 



0.018 - 

0.016 - 
0.014 

0.012 - 

I 

0.01 - 

0.008 

0.006 - 

0.004 - 

0.002 - 






without angular ordering Aqqq = 0.1 92 GeV 
with angular ordering Aqqq 0.221 GeV 
partons at hard collision 



0.2 



+ 
X 



0.4 



0.6 



■ I S i 

0.8 



Figure 7: The two x-spectra generated by the toy model, with angular ordering using 
Aqcd = 0.221 GeV and without angular ordering using Aqcd = 0.192 GeV respectively, 
look close to identical. The x-spectra are generated by backwards evolution from the scale 
of the hard collision, = m|, with a center of mass energy of 1800 Gev. For comparison, 
the x-spectrum for the partons at the hard collision is shown. 



where the same quark can be traced through the branchings back to the lower scale Q^. 
The same test as for forward evolution can be done, this time comparing only the x- 
distributions for valence quarks and doing the evolution from the higher scale to the 
lower Qq. Again the result is comforting. 

At this point the effect of angular ordering can be studied. The approximate opening 
angle of each emission is calculated using eq. (14) and an requirement, that the opening 
angle of an emission must be smaller than the previous one in the backward evolution, 
is introduced. This slows down the rate of emissions and consequently the rate of the 
evolution In x. To compensate for this Aqcd can be raised. To quantify the increase in 
Aqcd needed to compensate for angular ordering, x distributions is generated by back- 
wards evolution from an higher scale = m| (to be able to compare to Z° production 
in Pythia) down to a lower scale Q^. The backwards evolution is done with the center 
of mass energy set to 1800 Gev, as is the case for the Tevatron. At first the distribution 
is generated without the angular ordering requirement at the Aqcd value 0.192 Gev used 
in the tuning of the CTFQ5L parton densities, then distributions are generated with the 
angular ordering requirement at different values of Aqcd- The mean values of x can be 
plotted as a function of Aqcd and at some value of Aqcd the mean value of x is the 
same as without angular ordering using Aqcd = 0.192 Gev. According to the toy model 
algorithm, this new Aqcd needed to compensate for angular ordering is 0.221 Gev. As 
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Figure 8: A way of quantifying the increase in Aqcd needed to compensate for angular 
ordering for the process qq — > Z° at the Tevatron. First the x-spectrum for the initial 
partons is generated without angular ordering using Aqcd = 0.192 GeV. Then x-spectra 
for the initial partons are generated with angular ordering using different values of Aqcd • 
At Aqcd = 0.224 GeV the mean value of x is the same as without angular ordering at 
Aqcd =0.192 GeV. 



this value is obtained looking only at the mean value of x, it is necessary to check that 
the entire x-spectra look alike. This is done in Fig. 7. 

It should be noted that this study of the increase in Aqcd needed to compensate 
for angular ordering in the toy model only includes valence quarks in the backwards 
evolution. After this initial studies in the toy model, more thorough studies of the effects 
of the corrections to DGLAP evolution on the x evolution has been made using Pythia 
version 6.220, but based on the same principles. 

5.1.2 Pythia studies 

The same way of quantifying the increase in Aqcd needed to compensate for angular 
ordering can be used, with Pythia doing the backwards evolution. If this is done for Z° 
production at the Tevatron the plot in Fig. 8 is obtained. According to Pythia the Aqcd 
needed is 0.224 Gev, not very far from the value given by the toy model. As in Fig. 7 the 
shape of parton densities also agree well. 

This method can be used to find the Aqcd values needed to compensate for all three 
corrections for two different processes, both at the Tevatron and LHC. The Tevatron is a 
pp collider with a center of mass energy of 1800 Gev and the LHC is a pp collider under 
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construction with a planned center of mass energy of 14000 Gev. It is interesting to see 
how dependent the increase in Aqcd is on the center of mass energy and on the generated 
process. The processes, qq — > Z° and gg — > H°, are chosen to study the differences between 
having qq and gg in the initial state. H° denotes the Higgs boson, which is required within 
the Standard Model but has not yet been found. The Higgs mass is set to 120 Gev, which 
is somewhere in the neighbourhood of where it is expected to be. The results are presented 
in Table 1. 





Tevatron 


LHC 




qq^ Z° 


gg ^ H° 


qq^ Z° 


gg ^ H° 


angular ordering 


0.224 GeV 


0.239 GeV 


0.221 GeV 


0.236 GeV 


the condition that ii < 


0.26 GeV 


0.38 GeV 


0.21 GeV 


0.31 GeV 


as{p']_) instead of as{Q^) 


0.28 GeV 


0.30 GeV 


0.26 GeV 


0.27 GeV 


sum of all effects 


0.31 GeV 


0.50 GeV 


0.29 GeV 


0.46 GeV 



Table 1: The value of Aqcd needed to compensate for different corrections for qq — > Z° 
and gg — > H° at the Tevatron and LHC. 

The increase to compensate for angular ordering appears to be almost independent of 
the center of mass energy and not very dependent of which process is generated. It would 
of course be convenient to be able to use the same Aqcd value independently of the center 
of mass energy and particles in the initial state of the process studied. However for the 
other corrections the values seem to fluctuate more. 

Aqcd needs to be increased more when the initial state is gg, compared to having 
the initial state qq. This is not totally unexpected, since quark branchings and gluon 
branchings behave differently. The quark in the branching q — > qg tends to take a large 
momentum fraction. This, together with the ordering in Q^, means that the constraints 
on the DGLAP evolution are often fulfilled. In the branching g — > gg one of the gluons 
tends to take a large momentum fraction and the other a small momentum fraction. If 
the gluon of the continued evolution takes a small momentum fraction, the constraints 
are not as naturally fulfilled as for the case with quarks. It would therefore be possible 
to associate different effective Aqcd with different branching processes, but we did not 
pursue it here. There is also a difference in the increase in Aqcd needed between the 
different values of center of mass energy. This difference is not so easy to explain, but 
considerably smaller. 

5.2 The effect of an increased Aqcd value on p^-spectra and pri- 
mordial kj_ 

In Fig. 9 p_Lz-spectra for the process qq — > Z° at the Tevatron are shown. The p_Lz-spectra 
generated by Pythia, with Aqcd = 0.192 Gev and Aqcd = 0.31 Gev, are compared with 
experimental data [5]. The increase in Aqcd does not seem to affect the p_Lz-spectra much. 
To understand this it could be interesting to look at the effect of Aqcd on as at a typical 
scale of the parton shower, = 100 GeV^. When Aqcd is increased from 0.192 Gev 
to 0.31 Gev, as{Q^ = 100 GeV^) is increased by approximately 10%. The peak position 
of the p_Lz-spectra is shifted towards higher values by only approximately 1%. This is 
because the p_L-kicks caused by the emissions are in different directions in the xy-plane 
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Figure 9: p_Lz-spectra for the process qq — > Z° at the Tevatron, generated by Pythia 
and experimental data. The Pythia predictions are normalized to the total experimental 
cross section. 



(transverse to the direction of the beam), which means that they partly can cancel each 
other. If you instead look at the E± spectra, the sum of the absolute value of the p_L-kicks 
caused by emissions, the peak position is indeed shifted towards higher values by the 
expected 10%. 

Formally the p_L-kicks corresponds to emission of jets, so in theory the energy of these 
jets could be measured to get an experimental value for E± . Unfortunately, this alternative 
measure of activity in Z° events is less accessible experimentally, since also other physics 
could contribute to activity in an event. 

The absent effect of the increased Aqcd value on the p_Lz-spectra means that something 
else is needed to solve the problem with the too high value of the primordial k± needed 
in Pythia to fit experimental data. In Fig. 9 it is also obvious that in the p_Lz-spectra 
generated the event rates at the peaks are too high. In Subsection 5.3 another approach 
is described. 

For gg — > 11° at LHC the change in Aqcd has a bigger effect on the p_LH-spectra, see 
Fig. 10. However these spectra is not affected by primordial k± very much, which can be 
understood as follows. The center of mass energy at LHC is much greater, which means 
more emissions. At each branching the primordial k± is split between the two particles 
and not much is left at the production of the Higgs particle, cf. eq. (13). 
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Figure 10: pj_H-spectra for the process gg — > H° at LHC generated by Pythia. The Higgs 
mass is assumed to be 120 Gev. 



5.3 The effect of final-state showers of emitted particles with 
timelike virtuality on p^-spectra 

As seen in Subsection 4.4 final-state showers of emitted particles with timelike virtuality 
affect the p_L-spectrum of the produced particle. These timelike showers are supposed to 
have virtuality below the spacelike at the emission. It could be argued that this means 
they occur too late to influence the spacelike parton of the branching. The recoil would 
instead be taken by colour-connected neighbour partons. Since both Z° and H° are colour 
neutral, their p_L-spectra would then be unaffected. If this is right, the p_L-spectra can be 
simulated by switching off flnal-state showers of emitted particles with timelike virtuality. 
A study of the recoiling jets would obviously require more, but we leave that aside here, 
since we do not need to study these jets for the purpose of this thesis. In Fig. 11 this has 
been done and the agreement with experimental data is much better. 

For gg — > 11° at LHC the change in Aqcd has an even bigger effect on the p_L-spectra, 
when flnal-state showers of emitted particles with timelike virtuality are switched off. The 
peak position is then shifted towards higher vaiues by approximately 50%. 

In a new shower algorithm [13] the p_L-spectrum of the produced particle would not be 
affected by these timelike showers of emitted particles. This new algorithm is described 
in Subsection 5.4. 
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Figure 11: p_Lz-spectra for the process qq — > Z° at the Tevatron. Pythia generations with 
and without final-state showers of emitted particles with timelike virtuality, compared to 
experimental data. Both Pythia generations use Aqcd = 0.192 Gev and k± = 2 Gev. 



5.4 j9x-spectra generated by the new algorithm 

In the new algorithm the evolution is done m p\ = [1 — z)Q^ instead of Q^. This is 
possible since 

dpi (1 - z)dQ^ dg2 



pI {1-z)Q^ g2 (23) 

The ordering of emissions in p\ gives a set of colour connected partons. These partons can 
then evolve a final state shower as a system. The Z°-boson, as a colour neutral particle, 
is unaffected of colour dipoles. 

To compare evolution in p^ with evolution in Q^, the p_Lz-spectrum generated by the 
new algorithm can be compared to the p_Lz-spectrum generated by Pythia with final- 
state showers switched off. This is done in Fig. 12. The two spectra look similar except 
for small p±z where the new algorithm has better agreement with experimental data. 

To find which values of Aqcd and primordial k± that give the best fit to experimental 
data using the new algorithm, a test can be performed, is calculated by 



X 



2 _ (MC, - data,)2 

(72(MC,) + (T2(data,) 

where the summation is over the n different experimental data points. In Fig. 13 
values, generated with different values of Aqcd and primordial A;_l, are shown. The fit 



18 



30 



25 - 



20 



15 - 



10 



5 ^ 




experimental data i — \ — i 
Pythia, without fsr showers x 
Pythia, new algorithm 



X 



X 



-I 



X 



1^ 



10 

Ptz(GeV) 



15 



20 



Figure 12: The p_Lz-spectrum generated by the new algorithm compared to the p±z- 
spectrum generated by Pythia with final-state showers switched off. Both generations 
use Aqcd = 0.192 Gev and k± = 0.4 Gev. 



is best for low values of x^- A perfect fit would give a value equal to the number of 
degrees of freedom (the number of data points minus the number of parameters that are 
to be determined), which in this case is 48. Due to fairly large statistical fiuctuations and 
a lacking description of bin-to-bin correlations some of the values lie below this value. 
The values indicate that again the fit is best for Aqcd ~ 0.19 GeV and k± 2.0 GeV, 
but a trend can be seen that increasing Aqcd gives better fits for lower values of k±. 

In Fig. 14 p_Lz-spectra are generated with some Aqcd and primordial k± values that 
give good values. These p_Lz-spectra are generated using 10 000 000 events, instead of 
the 100 000 events used to generate the x^ values in Fig. 13. The x^ values then become 
46.5 for k^=2 Gev, Aqcd = 0.19 Gev, 51.4 for k^=lA Gev, Aqcd = 0.2 Gev and 74.6 for 
k± = 0.6 Gev, Aqcd = 0.22 Gev. In Fig. 14 it can also be seen that the two p_Lz-spectra 
generated with k± = 2 Gev, Aqcd = 0.19 Gev and k± = 1.4 Gev, Aqcd = 0.2 Gev 
give good fits to the experimental data, while the one generated with k± = 0.6 Gev, 
Aqcd = 0.22 Gev have too many events at low p±z and too few around the peak. 

Fig. 15 shows the same p_Lz-spectra as in Fig. 14, divided by the experimental data 
points, over a wider range in p±z- The p_Lz-spectra seem to differ mainly for small values 
of p_LZ, which is what is expected when different values of primordial k± are used. In 
Fig. 15, there is a trend that the p_Lz-spectra with lower values of primordial k± have 
more events with low p±z- A small effect of the change in Aqcd can also be seen for high 
values of p±z in Fig. 15, where higher values of Aqcd mean fewer events around the peak 
and more a bit above that. 
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Figure 13: foi" fits to experimental data using the new algorithm with different values 
of Aqcd and primordial A;_l 

At p_Lz values a bit above the peak many of the generated points he above the experi- 
mental data points. The fact that the points of the three different generations appear to 
fluctuate in unison owes to the point-to-point fluctuations in the data, where the statistics 
is much smaller than in the Monte Carlo studies. Again the two p_Lz-spectra generated 
with A;_L = 2 Gev, Aqcd = 0.19 Gev and A;_l = 1.4 Gev, Aqcd = 0.2 Gev flt the experi- 
mental data equally well, while the one with A;_l = 0.6 Gev, Aqcd = 0.22 Gev have too 
many events with low p\_z- However it should be noted that the points where the flt of 
the p_Lz-spectrum with A;_l = 0.6 Gev, Aqcd = 0.22 Gev is considerably worse than the 
others are mainly the three flrst. A systematic error in the experimental data could affect 
closely lying points in the same way and thus bias the conclusions. 

The new algorithm does not yet include heavy quarks (b, c). When p_L-spectra gener- 
ated by the ordinary Pythia algorithm are studied, the exclusion of heavy quarks does 
not have any measurable effect, however. 

6 Inclusion of heavy quarks in the new algorithm 

This section is an attempt to flnd a way of introducing heavy quarks in the new algorithm. 
No obvious way to do this is found and the work is under development. This is only a 
beginning of these studies. 

6.1 q ^ qg 

To derive the probability for a quark to emit a gluon, the cross section for a process 
that involves a quark emitting a gluon can be divided by the cross section for the same 
process without the emission of the gluon. The process used in this derivation, shown in 
Fig. 16, involves a made-up Higgs particle that couples to the particles needed to make 
the processes as simple as possible with a vertex factor of 1. This is allowed since the aim 
is to derive the probability of the quark emitting a gluon in the colinear limit where it 
does not depend on the rest of the process. 

The needed matrix elements squared are calculated using CompHep [14]. One prob- 
lem of technical character, is that CompHep require the Higgs to be unstable, so we 
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Figure 14: p_Lz-spectra generated by the new algorithm with Aqcd and primordial k± 
values that give good values. 



have to generate an isotropic Higgs two-body decay. This makes the calculations more 
complicated. 

6.1.1 Derivation of the DGLAP splitting function for massless quarks 

The cross section for the process in Fig. 16b is given by 



(7 



dxidx2f{xi)f{x2) ^sdPSs (25) 



where Xi, X2 are the fractions of the total proton momentum carried by the initial quarks, 
/ is the structure function of the proton (neglecting the dependence), dPSs is the 
integration over the phase space for the three final particles and the constituent cross 
section is given by 



|2 



|M2^3| r/ 2 ■2\ /r,r.\ 

0-3 oc — d{m- - m^j (26) 



where |M2_>3|^ is the matrix element squared and mi is the mass of the intermediate 
particle. 

With the particles denoted according to Fig. 17, the phase space integral is given by 
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Figure 15: The same p_Lz-spectra as in Fig. 14, divided by the experimental data points 
and over a wider range in p±z- The experimental errorbars are also divided by the 
experimental data points. 



where P denotes four-momentum. This can be rewritten as 

dPSs = dPS2 • — ^ / dmf ^ (28) 

4 • [ZTT)'' J X1X2S 

where dPS2 is the phase space integral for the process in Fig. 16a and t is the Mandelstam 
variable for the process qq — > H°g. The cross section given in eq. (25) then becomes 



a oc 



1 f |Af"2->3|^ dt 

A TTT^ h^^i^^^fi^i)fi^^) ^r?^dmf dPS2 5(mf-mH) (29) 

4 • [Ztt j J Zs X1X2S 

In the parton shower approximation, the cross section for the process in Fig. 16b can 
also be calculated using the cross section for the process in Fig. 16a and the DGLAP 
splitting function for q — > qg 



croc /da;ida;2/(xi)/(a;2)-^^%i^d2;^^^Pq^qg(2;)dPS2^(xia;2^;s - m^) (30) 
J Zs Lj Z'K 

In the coUinear limit where 

m\ = X1X2ZS =^ dm? = XiX2sdz (31) 

eq. (29) becomes 
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H° 





H° 




Figure 16: a) and b) qq — > H° — > e with and without the emission of a gluon. H° is 
a made-up scalar particle that has the couplings necessary for the processes. 




Figure 17: The final particles of the process in Fig. 16b. 3 is the gluon, i is the intermediate 
made up Higgs boson, 1 and 2 are the e"'"e~-pair and denotes the available energy or 
momenta from the collision. 



1 



G OC 



4 • (2^1 



\M< 



2->3 



2s 



■ dz dt dPS2 ^{xxx^zs — m^) 



(32) 



The condition that the two expressions for the cross section, eq. (30) and eq. (32), should 
be equal gives, after some further considering 



q->qg 



(-) 



2->3 



(33) 



The matrix elements squared can be calculated, e.g. using CompHep. The process in 
Fig. 16b can not be separated from the process where the other quark emits the gluon. 
The matrix elements must be added and then squared. In the coUinear limit where the 
gluon is emitted in the direction of the first quark the matrix element squared of the 
second process disappears. The interference term however can contribute to the derived 
expression for the splitting function. Using the CompHep matrix elements squared gives 
the following expression for the splitting function 



q->qg 



(-) 



4 ^2(^2^ 11 



3s(s(l -^) + 



(34) 



In the coUinear limit, where is small compared to s, eq. (34) goes to the DGLAP 
splitting function in eq. (2). 
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6.1.2 Massive quarks 



The corresponding calculation can be done for the case when the quark emitting the gluon 
in the process in Fig. 16b has a non-negligible mass. For a c quark the calculation gives 



(3^2^2 ^ ^2Q2 _ 2znil + m2 + Q'')s - 2mlQ'' - 2mlz 

(35) 



dPc^cg a, 4 {?,z^ml + z^Q^ - 2zml + + Q^)s - 2mlQ^ - 2mlz 



dg2 2^3 {ml + Q^f{{l 

Since the new algorithm is ordered in what is needed is really ^'"^.^^^ instead of 
For the case of a c quark emitting a gluon function of becomes 





ml){l-z)-Q^)Q^s 




[s — miy 


Q' = 


{s-ml){l-z) 




2 




1 s — ml 







(36) 



Neglecting the Q"^ term gives 



m1)p\ 



{l-z)s 
dQ"^ s — ml 



dp\ (1 — z)s 



(37) 



Fven with these simplified expressions for g2 and the expression for gets very 

complicated. 

Complexity in itself is not an absolute hindrance, however: if a simple upper estimate 
can be found, standard rejection techniques can be used to obtain results in agreement 
with the full expression. It therefore remains to find some suitable such estimate in a 
p_L-oriented framework. 



6.2 g ^ qq 

To derive the probability for a gluon to split into a cc pair, the same method as in 
Subsection 6.1 can be used, this time with the processes shown in Fig. 18. The following 
expression is obtained when, again using CompHep, the matrix element squared for the 
process in Fig. 18b is divided by the matrix element squared for the process in Fig. 18a 

(2z'mg + 2z'Q' - 2zQ'' + mg + Q'')s 2zQ^ - 2mgQ' - 2Q^ - 2m^ - 2m^z 

^ (ml + Q'^Y is — ml) ^ (ml + Q'^V (s — mD 



m 



slmf: + (J2)2 ig — fji2\ 
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H° 





H° 




Figure 18: qq — > H° — > e . In b) the quark comes from a gluon splitting into a qq-pair. 
is again a made-up scalar particle that has the coupling necessary for the processes. 



If the fact that there is a massive c quark in the initial state of the process in Fig. 18a 



is taken into account the denominator in Fq. 26 should be 2(s 



instead of 2s. This 



means that the expression in Fq. 38 should be multiplied with ^^-4^, which gives the 
following expression 



47ra« 



+ (1 + 2^2)^2 



{ml + g2)2 



s(m2 + g2)2 + 



(39) 



where the first term takes the form of the DGLAP splitting function when m\ is set to 0. 

In the process in Fig. 18b there is a massive c quark in the final state. This means 
that the phase space integral gets more complicated 



dPSa = dPS, 



1 



4 • (2^1 



\J{xiX2S — mf + m^)2 — AxiX2sm1 
X1X2S — m? + ml 



dm? 



dt 



X1X2S 



(40) 



Again it remains to find a simple upper estimate as a starting point for a more precise 
inclusion of mass effects. 



7 Summary 

To explore perturbative and nonperturbative effects in transverse momentum generation 
p_L-spectra have been studied. More specificaly, the infiuence of corrections in Pythia to 
pure leading log DGLAP evolution on x evolution and generation of pj_-spectra has been 
investigated. This has been done in hope that the reason for the need of a primordial A;_l 
much larger than expected might be found. 

To compensate for the corrections to pure leading log DGLAP evolution in Pythia, 
Aqcd needs to be raised to the values given in Table 1. For qq — > Z° at the Tevatron 
this increase in Aqcd does not affect the generated p_Lz-spectra very much. This can be 
understood as the effect of p_L kicks of the branchings in different directions in the plane 
transverse to the beam direction partly canceling each other. The generated p_Lz-spectra 
do not fit the experimental data very well. For gg — > 11° at LHC a changed Aqcd has 
a larger effect on the generated p_LH-spectra, due to higher center of mass energy and 
consequently more branchings. In this case the primordial A;_l of the initial partons does 
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not affect the generated p_LH-spectra very much, since it is spiit between the partons of 
the many branchings. 

Final-state showers of emitted particles with timelike virtuality affect the generated 
p_L-spectra. It could be argued that these showers occur too late to influence the spacelike 
parton of the branching, since they are supposed to have virtuality below the virtuality 
of the branching. When these are turned off the fit to the experimental data for qq — > Z° 
at the Tevatron is much better. 

In a new shower algorithm the branchings are ordered in p\ instead of Q^. This gives 
a set of colour connected partons which can then undergo final state showers as a system. 
These final state showers would not affect the produced Z°, since it is a colour neutral 
particle. In this algorithm the mcrease m Aqcd has a larger effect on the generated 
p_LZ-spectra. When three different p_Lz-spectra, using values of Aqcd and primordial k± 
that give good values, are studied the fit is equally good for the two spectra using 
k± = 2 Gev, Aqcd = 0.19 Gev and k± = lA Gev, Aqcd = 0.2 Gev respectively, while the 
fit is a bit, but not very much, worse for the one using k± = 0.6 Gev, Aqcd = 0.22 Gev. 
The problem with a large value of primordial k± is not completely solved but reduced. 

The studies to include heavy quarks in the new algorithm are under development. No 
obvious way of doing this has yet been found. 
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